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^ ■ Abstract. We prove that the 3-D compressible Euler equations with surface tension 

along the moving free-boundary are well-posed. Specifically, we consider isentropic dy- 
namics and consider an equation of state, modeling a liquid, given by Courant and 
Fricdrichs [8] as p(p) = ap 1 — ft for consants 7 > 1 and a, (3 > 0. The analysis is made 
difficult by two competing nonlinearities associated with the potential energy: compres- 
sion in the bulk, and surface area dynamics on the free-boundary. Unlike the analysis of 
the incompressible Euler equations, wherein boundary regularity controls regularity in 
the interior, the compressible Euler equation require the additional analysis of nonlinear 
wave equations generating sound waves. An existence theory is developed by a specially 
chosen parabolic regularization together with the vanishing viscosity method. The ar- 
tificial parabolic term is chosen so as to be asymptotically consistent with the Euler 
equations in the limit of zero viscosity. Having solutions for the positive surface tension 
problem, we proceed to obtain a priori estimates which are independent of the surface 
tension parameter. This requires choosing initial data which satisfy the Taylor sign con- 
dition. By passing to the limit of zero surface tension, we prove the well-posedness of the 
compressible Euler system without surface on the free-boundary, and without derivative 
loss. 
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1. Introduction 

1.1. The compressible Euler equations in Eulerian variables. The compressible Eu- 
ler equations with moving free-boundary are given by the following system: 



d t {pu) +div(jOW(giu+pId) = in Q(t), (1.1a) 

d t p + div(pu) = in fi(t), (1.1b) 

p = aH(t) on (1.1c) 

V(T(t))=u-n(t) (l.ld) 

(u,p) = (u ,p ) on O(0), (Lie) 

n(o) = n, (l.if) 



where £l(t) denotes an open and bounded subset of K 3 , T(t) = dfl(t) is the moving free- 
boundary, and t € [0,T] denotes time. We use the notation V(T(t)) for the normal velocity 
of boundary T(t), which is equal to the normal component of the fluid velocity u ■ n, where 
n(t) is the outward-pointing unit normal to T(t), u — (u\, U2,u^) denotes the velocity field, 
p denotes the pressure, and p denotes the density. 

The first two equation are conservation laws for momentum and mass. The boundary 
condition (|l.lc[) is often referred to as the Laplace- Young condition, stating that the fluid 
stress is proportional to the mean curvature H(t) of the moving surface, the proportionality 
constant defining the surface tension parameter a. The last two equations provide the initial 
conditions for the dynamics. 

In order to model the motion of a compressible liquid, we use the equation-of-state given 
by Courant and Friedrichs [5] as 

p(x, t) = ap(x, t) 1 - /3 for 7 > 1, (1.2) 

where a > and /3 > 0. For convenience, we set a = 1. 

Using the equation of state (|1.2I) . the momentum equations (|l.la[) and Laplace- Young 
boundary condition (jl.lcl) are equivalently written as 

p[d t u + (u ■ D)u] + Dp~< = infi(t), (1.3a) 

pi=f3 + (jH onT(t). (1.3b) 

We assume that the initial density function is strictly positive and that 

Po > A > in Q. 

In the absence of surface tension, we further require the initial pressure function p$ to satisfy 
the Taylor sign condition (see, for example, [30] and [26]), given by 

where N denotes the outward unit normal to T. This is equivalent to 

Fin 1 

< v < -^P- on T. (1.4) 

on 

1.2. Prior results on the Euler equations with moving free-boundary. 



^Using 11.21 . liquid water is modeled using the values 7 = 7, a = 3001 and /3 = 3000. 
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1.2.1. The incompressible setting. There has been a recent explosion of interest in the anal- 
ysis of the free-boundary incompressible Euler equations, particularly in irrotational form, 
that has produced a number of different methodologies for obtain- ing a priori estimates. 
The accompanying existence theories have relied mostly on the NashMoser iteration to deal 
with derivative loss in linearized equations when arbitrary domains are considered, or on 
complex analysis tools for the irrotational problem with infinite depth. We refer the reader 
to [^ITT | [T4 l [2T 1 [2 ^ [27 l l35 l [36 l l39] for a partial list of papers on this topic. 



1.2.2. The compressible setting. The mathematical analysis of moving hypersurfaces in the 
multidimensional compressible Euler equations began with the existence and stability of the 
shock- front solution initiated in [24] and extensively studied by [16l - [T8|[2"5] (see the references 
in these articles for a thorough review of the literature in this area.) More delicate than 
the non-characteristic case of the shock-front solution, the characteristic boundary case is 
encountered in the study of vortex sheets or current vortex sheets. This class of problems 
has been studied by [¥,6,7,32,33] (and see the references therein); the linearization of the 
vortex-sheet problem produces derivative loss, similar to that experienced by many authors 
in the incompressible flow setting (both irrotational flows and flows with vorticity). 

The problem of the expansion of a compressible gas with the so-called physical vacuum 
singularity has been studied in j9[[T^[T3][T9l[20] , and is degenerate because of the vanishing 
of the density function on the moving free-boundary. 

For the model of a compressible liquid, considered in this paper, the Euler equations 
are uniformly hyperbolic thanks to the equation-of-state (|1.2[) . In the absence of surface 
tension, an existence theory was given in |23j using Lagrangian coordinates and a Nash- 
Moser construction, but the estimates had derivative loss. Using the theory of symmetric 
hyperbolic systems, the paper [33] gave a different proof for the existence of solutions, also 
with derivative loss. We prove well-posedness for the motion of a compressible liquid with 
and without surface tension, and with no derivative loss. We also establish asymptotic limit 
of zero surface tension. 

1.3. Fixing the domain and Lagrangian variables. To transform the system (|1.1[) into 
Lagrangian variables, we let r)(x, t) denote the flow of a fluid particle x at time t. Thus, 

dtf] = u o i] for t > and T)(x, 0) = x 

where o denotes composition, so that [u o r]](x,t) = u(j](x,t),t). The flow map rj induces 
the following Lagrangian variables on Q: 

v = uo rj (Lagrangian velocity), 

f = p o rj (Lagrangian density) , 

A = [Dn]" 1 (inverse of the deformation tensor), 

J = Act Dn (Jacobian determinant), 

a = J A (cof actor matrix of the deformation tensor). 

Using the notation defined below in Section l2.1.2l the Lagrangian version of equations (|1.1[) 
on the fixed domain Q is given by 

fvi + A>?r, k = inOx(0,T], (1.5a) 

ft ■ /,!;<''., inOx(0,T], (1.5b) 

P=P + aH(rj) onrx(0,T], (1.5c) 

(v> v >f)\t=o = (e,u ,p ) onfi, (1.5d) 
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where T — dtt and e(x) = x denotes the identity map on fi. 

Since J t — a\v l ,j by (|2.5j) . we multiply (|1.5bj) by J and integrate in time for the identity 

f = PoJ- 1 - (1.6) 

For a > and 7 = 2, we use the identity (|1.6[) for the Lagrangian density / and equivalently 
write the compressible Euler equations (| 1 . 5|) as 

Pov] + a\ (p 2 J" 2 ), fc = in fi x (0, T], (1.7a) 

pp~ 2 = p + aH{r]) onrx(0,T], (1.7b) 

(rj,v)\ t =o = (e,u ) on Q. (l-7c) 

For <7 > 0, we shall refer to the equations (11.71) as the surface tension problem. 
For reference, we record that the momentum equations (|1.7aj) are equivalent to 

v\ + 2A k i {p Q J- 1 ), k =Q inQx(0,T]. (1.8) 

Remark 1 . The identity (jl.6l) establishes that in Lagrangian variables the initial density po 
is a parameter in the system of compressible Euler equations. 

1.4. The higher-order energy functions E(t) and £(£). While the physical energy 
Jnl^oiM 2 + PqJ' 1 + @J\ + o~A(t), A{t) denoting the surface area of T(t), is a conserved 
quantity, it is far too weak for the energy estimates methodology that we employ. We instead 
define the higher-order energy functions E(t) and £(i) to respectively correspond with the 
surface tension problem and the zero surface tension limit. Although neither E(t) or £(i) 
is conserved, we will establish that each of sup tg [ T j E(t) and sup tS [ T i £(t) is bounded on 
a sufficiently small time-interval of existence [0,T]. 

1.4.1. The higher-order energy function for a > 0. We define the energy function E(t) as 

5 2 

E(t) = 1 + J2 WdtV(t)\\l-a + \vttt ■ n{t)\\ + J2 \d 2 d?v ■ n(t)\i 5 _ a . (1.9) 

o=0 a=0 

We let Mq > denote a generic constant given by a polynomial function P of E(0): 

M o = P(E(0)). (1.10) 

1.4.2. TTie higher-order energy function for a = 0. We define the energy function £(i) as 

£(<) = l + J2 W d Mt)\\l. 5 - ia + E 11^ J Wll4.5-Ia + II (1-H) 
a=0 a=0 

We let Mo > denote a generic constant given by a polynomial function P of £(0): 

M o = P(£(0)). (1.12) 

Section [2~T1 explains the notation in (|1.9I) and (jl.ll|) . 

Remark 2. Corresponding with cr = 0, time-derivatives in the higher-order energy function 
£(i) scale like one-half a space-derivative. Since the scalar product of the momentum equa- 
tions (|1.7ap and the tangential vector r), a yields the identity poVt ■ i], a = —J{p 2 ] J^ 2 ), a in O, 
the Laplace- Young boundary condition (jl.7bl) with a — provides that 

vt ■r), a = on T. 

Differentiating this identity with respect to time shows that dt scales like {d x ) x / 2 in the zero 
surface tension limit of (11.71). 
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1.5. The initial data (po,uo,il). We assume that po and uo are given and sufficiently 
smooth. We require that the initial density po satisfy 

p >2A>0 inO. (1.13) 

Using the identities 77(0) = e and J t = a s r v r , s , we let J\ be defined by 

Ji = <9 t (J~ 2 )(0) = -2divu . 

We let Vi and vn be the vectors respectively given by 

v{ = -2d iPo and vi = -po 1 [d i {p 2 J 1 )+d t a'?{0)(p 2 ), k ], 

where dta(0) is a smooth function of Duq. We define, as a function of po and uq, 

J a =d?{J- 2 )\ t =o for a = 0,1, 2, 3. 

1.5.1. The case of positive surface tension. For a > 0, we assume that Q is an H -class 
domain and that po and uq are in H 4 (Q). We define, as a function of po and Uo, 

H a =d?H{r,)\ t=0 for a = 0,1, 2, 3. 

We require that the initial data satisfy the following compatibility conditions: 

p 2 J a = d?/3 + aH a on V for a = 0,1,2,3. (1.14) 

For [3 as in (|1.2p and A as in (jl . 13|) , we note that a = in (|1.14[) yields 

crff > 4A 2 - (3. 

1.5.2. The case of zero surface tension. For a — 0, we assume that f2 is an _ff 4 5 -class domain 
and that po is in £f 4 - 5 (fi) and u is in iJ 4 (51). 

We require that po and ^ satisfy the Taylor sign condition: 

0<2v<~^Walat{plJ- 2 ), k \ t=0 onT. (1.15) 

We require that the initial data satisfy the following compatibility conditions: 

p 2 J a = d?/3 on T for a = 0,...,6. (1.16) 

Remark 3. The twice-mean-curvature function H (77) is not present in the zero surface tension 
limit of (jl.7l) . Accordingly, there is no restriction on the curvature of the initial surface T. 

1.6. Main Results. The main results of this paper are the existence and uniqueness of 
solutions to the surface tension problem and its zero surface tension limit. 

Theorem 1.1 (Existence and uniqueness for a > 0). Suppose that the initial data (po, uq, £1) 
verify 

(1) M a = P(E(0)) < 00, 

(2) the lower-bound condition (| 1 . 1 3|) . and 

(3) the compatibility conditions (|1.14p . 

Then for some T > 0, there exists a solution to (|1.7[) on the time-interval [0,T] such that 
p(t) > A in Q(t), aH(t) > -/3 on T(t), and 

sup E(t) < 2M . 
te[o,T] 

Furthermore, the solution is unique if the initial data is such that 

E W d ?n(0M-a + hut ■ n(0)\l + £ \d 2 d?v ■ n(0)| 2 . 5 _ Q < 00. 
a=0 a=0 
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Theorem 1.2 (The zero surface tension limit). Suppose that the initial data (po,uo,Q) 
verify 

(1) M = P(£(0)) < oo, 

(2) the lower-bound condition (|1.13l) , 

(3) the Taylor sign condition (|1.15[) . and 

(4) the compatibility conditions (|1.16[) . 

Then for some T > 0, there exists a solution to (|1.7p tuii/i cr = on the time-interval [0, T] 
suc/i i/iai p(i) > A in Q(t), — > v on T(t), and 

sup E(t) < 2M . 
te[o,T] 

Furthermore, the solution is unique if the initial data is such that 

E Wd^M.^a + E ll9M0)llL-| O + ll^(o)iii < oo. 

o=0 a=0 

Remark 4. The proofs of Theorems 11.11 and 11.21 do not rely on our choice of 7 = 2 and as 
such are valid for general 7 > 1, since by introducing new enthalpy- type variables for p^~~ 
and J 7-1 in (|1.5[) . we can reduce the case of general 7 > 1 to the case that 7 = 2. 

1.7. Structure of the proofs of Theorems 11.11 and 11.21 

1.7.1. Existence for the surface tension problem (|1.7[) . An existence theory is obtained via 
the vanishing viscosity method, but with a very special choice of artificial viscosity which 
does not alter the transport-type structure of vorticity. This is introduced in Section [3] The 
Euler equations (|1.7a|) yield 

curLj Vt = 0, 

where the Lagrangian curl operator curl I( is defined below in Section 12.1.31 By defining 
a parabolic approximation of the Euler equations (|1.7a|) which preserves the homogeneous 
vorticity equation, we ensure that the vorticity estimates for the approximate problem are 
independent of the parabolic approximation parameter re. The parabolic approximate re- 
problem defined in Section [3] maintains a homogeneous vorticity equation. 

The structure of the Euler equation (ll.7ap shows that an a priori estimate for vt provides 
an estimate for the gradient of J, which then yields an estimate for the gradient of divry. 
As such, yet another constraint on the choice of parabolic regularization operator is that its 
presence must still permit estimates for D div rj whenever estimates for vt are known. We 
choose our re-parabolic operator so that the parabolically-regularized momentum equations 
have the form 

/ + «/* = 9, 

for / equalling the gradient of J. Such a structure allows us to obtain re-independent 
estimates for /, given estimates for g. We must further add parabolic regularization to the 
surface tension boundary conditions (|1.7bj) which also has the same structure, in order to 
infer the maximal regularity of the free-boundary. 

Our strategy is to obtain energy estimates for the highest number of time-derivatives in 
the problem, and then use the structure of the momentum equations to infer estimates for 
the divergence. In conjunction with curl-estimates and boundary regularity, we boot-strap 
the full regularity of the velocity, its time-derivatives, and the flow map rj. 

In Section 21 we perform energy estimates in the fourth (highest) time-differentiated 
approximate re-problem. From these energy estimates, we are able to close estimates for 
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vtttt(t) in L 2 (VL) and vttt{t) in We then infer the optimal regularity of vtt and vt- 

This, in turn, enables us to get the optimal regularity for v and r\. All estimates arc found 
to be independent of the artificial viscosity parameter n. 

In Section [5l we construct smooth solutions the approximate k- problem. We utilize (the 
Lagrangian variables version of) the basic vector identity — A = curl curl— Ddiv in order to 
replace gradient of the divergence of v with the Laplacian of v, modulo lower-order terms. 
Additionally, we use that a sufficiently regular vector £ € M 3 satisfies 

^'4A^,fe=^J- 1 (curl n x n + ^J- 1 (div„£)n + &(£,7 7 ), 

where &(£, 77) is a first-order differential boundary operator with respect to £ and a second- 
order differential boundary operator with respect to r\. Since the regularity of the flow map 
77 is dictated by the regularity of v in our construction scheme, we employ a horizontal 
convolution-by-layers operator (first introduced in in order to view b(v,r]) as a lower- 
order term; with the horizontal convolution-by-layers approximation in place, solutions can 
be found via the existence theory for uniformly parabolic second-order equations. We must 
then find estimates for such solutions which are indeed independent of the convolution 
parameter and pass to the limit. 

In Section [51 we use the K-independent a priori estimates established in Section @] and 
the construction of solutions to our approximate ^-problem in Section [5] to establish the 
existence and uniqueness of solutions to the surface tension problem (jl.7[) . 

1.7.2. The zero surface tension limit. In Section [71 we establish our existence theory for 
the zero surface tension limit of (|1.7[) via u-independent a priori estimates. For initial data 
satisfying the Taylor sign condition (|1.15|) . we have that solutions to the surface tension 
problem (|1.7|) satisfy 

< \V ■ n\ 2 < jf ^4a k MJ-%k \V • n\\ 

on a sufficiently small time-interval [0, T]. 

2. Preliminaries 

2.1. Notation. 

2.1.1. The three-dimensional gradient vector. Throughout this paper the symbol D will be 
used to denote the three-dimensional gradient vector 

D=( —,—,—1 . 

\dx\ 8x2 dx3 J 

2.1.2. Notation for partial differentiation and Einstein's summation convention. The fcth 
partial derivative of F will be denoted by F,k = J^-. Repeated Latin indices i,j, k, etc., are 
summed from 1 to 3, and repeated Greek indices a, /3, 7, etc., are summed from 1 to 2. For 
example, F nl = £? =1 and F\ a I^G\ = ELi E? =1 f£^ Q/3 g- 

2.1.3. The divergence and curl operators. We use the notation div u to denote the divergence 
of a vector field u on fl: 

divu = u ,1 +m 2 ,2 +u 3 ,3 ) 
and we use the notation curl it to denote the curl of a vector u on Q: 
curl it = (u 3 , 2 -w 2 ,3 ," 1 ,3 -m 3 ,i ," 2 ,i -u 1 ^) ■ 
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We recall that the permutation symbol Eijk is defined by 

!1, even permutation of {1, 2, 3}, 
— 1, odd permutation of {1, 2, 3}, 
0, otherwise. 

This allows for the curl of a vector field u to be expressed as curlu = e.jkU k ,j. Letting 
v = u(r/) for a given flow map rj, we use the notation div,, v to denote the Lagrangian 
divergence of v on f2: 

div,, v = A s r v r , s , 

and we use the notation curb, v to denote the Lagrangian curl of v on fi: 

curlew = e. 3k AjV k , s . 

2.1.4. The scalar- and cross-product of vectors in R 3 . Let u and v be vectors in R 3 . The 
scalar-product of u and v, denoted u ■ v, is defined as 

u ■ v = uV + u 2 v 2 + u 3 v 3 . (2.1) 

The cross-product of u and v, denoted u x v, is defined as 

u x v — e. 3 ku j v k . (2.2) 

2.1.5. Local coordinates near T. We let Q C R 3 denote an open, bounded subset of class 
H s for s > 4, and we let {Ui]f =1 denote an open covering of V — dil, such that for each 
I G {1,2,...,K}, with 

Bi = B(0, ri), denoting the open ball of radius n centered at the origin, 

B+ = Bi n {x 3 > 0}, 

T>i = Bif) {x 3 = 0}, 

there exists an iJ s -class chart Q\ satisfying 

9i : Bi —> Ui is an H s diffeomorphism, 

0i(Bf) = Uinn, 
e l (D l ) = u l nv. 

For L > K, we let {Ui}f =K+1 denote a family of open balls of radius r\ properly contained 

[ UkA }\ 

Uk Uk+i 

...y ih 



Figure 1. Indexing convention for the open cover {Ui}[ =1 of 51. 
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in Q, such that {Ui}f =1 is an open cover of il. We let 

denote a C°° partition-of- unity subordinate to the open covering of f2. 

2.1.6. Tangential derivatives. On each Ui n fi, for 1 < I < L, we let Bi denote the local 
tangential-derivative. That is, for a differentiable function / on ft, the ath component of 
the local tangential-derivative of / is defined in Ui n by 

W=(^[/oft])o«r 1 = (woft)^)o-?r 1 , 

where for K + 1 < I < L, we have set 0/ to be the identity map e. 

We let 8 denote the tangential-derivative whose ath component is given by 

L 
i=l 

We use <9 Q / or /, a to denote the components of the tangential-derivative of /. 

2.1.7. Geometry of the moving surface T(t). The vectors n, a for a = 1, 2, span the tangent 
space to the moving surface = r)(T) in R 3 . The surface metric g on r(i) has components 

g a ji = Via -ViP ■ 

We let go denote the surface metric of the initial surface T. The components of the inverse 
metric [g]^ 1 are denoted by g a " . We use ^fg to denote ^ det g; we note that y/g = xry,2 |, 
so that n = [n,i xn,2]/^/g. Equivalently, 

Vgn = n, 1 xn, 2 - (2.3) 
The Laplace-Beltrami operator A g is dehned on T as 

2.1.8. Sobolev spaces on il. For integers k > and a smooth, open domain 51 of R 3 , we 
define the Sobolev space H k (fl) (# fc (0;R 3 )) to be the closure of C°°(n) (C°° (H;M 3 )) in 
the norm 




for a multi-index a e , with the convention that \a\ = a\ + a 2 + a s . For real numbers 
s > 0, the Sobolev spaces H s (tt) and the norms || • || s are defined by interpolation. We will 
write fP(fi) instead of _ff s (0; R 3 ) for vector- valued functions. We use H s (fiy to denote the 
dual space of H s (£l). 

2.1.9. Sobolev spaces on T. For functions u € H k (T), k > 0, we set 




for a multi-index a G For real s > 0, the Hilbert space H S (T) and the boundary norm 
| • | s is defined by interpolation. The negative-order Sobolev spaces H~ S (T) arc defined via 
duality. That is, for real s > 0, 

H- S (T) = H S (T)'. 
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Remark 5. Throughout this paper, we suppress the Euclidean measure dx by letting J n 
represent J n dx. Similarly, the notation J r represents J r dSa where dSo = ^fg^dxxdx^ is 
the surface measure of the initial surface T. Equally, the time integral J* should be read as 

Jo ds - 

Remark 6. We let | • \ s .t> 1 denote the if s (2?;)-norm. 

2.2. Differentiation and geometric identities and properties. 

2.2.1. An identity for the Jacobian determinant J. With dim 51 = 3, we have that the 
Jacobian determinant J is written as 

J= dimn^]' W 
which follows from the definition of the cofactor matrix a in Section 11.31 



2.2.2. Time- differentiating the Jacobian determinant J and the cofactor matrix a. We record 
the following basic differentiation formulas: 

d t J = a s r v r , s , (2.5) 
dta'y = J- 1 [a s r a'y-a!a k r }v r , s . (2.6) 

Using (|2.5p and (|2.6p and the fact that a — J A, we have that 

d t Al = -A\A k r v r , s . (2.7) 

We note that the time-differentiation formulas (|2.5p ~ (|2.7p at once become formulas for 
tangential-differentiation by replacing v r , s with dn r , s in the right-hand sides of (|2.5p ~ (|2.7p . 
The formulas (|2.5p and (|2.6p imply the following scaling relations: 



dt J ~ aDv, dta ~ a 2 Dv, 

d 2 t J - a 2 (Dv) 2 + aDv t , d 2 t a - a 3 (Dv) 2 + a 2 Dv t . 

These scaling relations are particularly useful when estimating error-terms. 

2.2.3. The Piola identity. Columns of every cofactor matrix are divergence-free. Thus, 

atk = 0- (2.8) 

2.2.4. Relating the normal vectors of T and T(t). With TV as the outward unit normal to 
the reference surface T, the outward- normal direction of the moving surface T(t) is 

a k t N k = |77,! X77, 2 |n\ 



The identity y/g — ydet g = \r],i x?7, 2 | implies that 

a k N k = ^gn\ (2.9) 
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2.2.5. Derivatives of the inverse metric g a " , Jacobian ^fg and unit normal n. A tangential- 
derivative of the inverse metric g al3 , Jacobian determinant ^[g and moving outward normal 
unit vector n are given by the formulas 



Also, 



B g <*0 = -g^dg^g"?, 


(2.10a) 


<V# = ^y/gg^dg^, 


(2.10b) 


Bn = -g~t s [8ri, s -n]?7, 7 . 


(2.10c) 


dtg afi = -g^dtg^, 


(2.11a) 


dty/g = ^y/gg^dtg^v, 


(2.11b) 


d t n = -g jS [v, s "n]r) n . 


(2.11c) 



Remark 7. The right-hand side of ()2.10cj) and (|2.11cj) is a vector that is tangent to the 
embedded surface. 

2.2.6. Relating the Laplace- Beltrami operator A g to the unit normal n. With the formulas 
(|2.10al) and (|2.10b[) . we have that the Laplace-Beltrami operator A g = y/g^dal^fgg^ Bp] 
applied to the particle flow 77 decomposes into normal and tangential components as 

g af3 bhai3-n]n 

We therefore have the identity 

y/gA B (<r,) = V99 aP h a -n]n. (2.12) 
For reference, we recall the identity A g (ri) = —H(rf)n. 

2.3. Two identities for the Euler equations in Lagrangian variables. 

2.3.1. The Lagrangian vorticity equation. With the operator curb, defined in Section 12.1.31 

curl^ Vt — in Q. (2.13) 
The identity (|2 . 13[) is obtained by taking the Lagrangian curl of the Euler equations (11.81) . 

2.3.2. A tangential identity for vt on the boundary T. Setting a = 1, 

VfV,a=f~ 1 d a (g flv [n,^-n]) onT. (2.14) 

The identity (|2.14l) is established using the Euler equations (|1.8I) , the Laplace- Young bound- 
ary condition (|1.7b|) . and the formula (I2.12j) . 

2.4. General inequalities. 

2.4.1. Trace estimates. For s > i and some constant C independent of w € H s (il), the 
trace theorem [TJ states that the trace of w is defined in H S ~^(T) with the estimate 

\w\ s _ h <C\\w\\ s . 

Lemma 2.1. Let w e L 2 (0,T; iJ 1 (il)) n L°°(0, T; L 2 {Q,)). Then 

I Ml M <&l \\M\x + C 5 T sup \\w{t)\\l (2.15) 

JO Jo t£[0,T] 

where the constant C$ depends on 1/5 > 0. 



12 



D. COUTAND, J. HOLE, AND S. SHKOLLER 



Proof. Using interpolation and Young's inequality with S > 0, we have that 

\\w\\l 75 <C [ T \\w\\l\\w\\f <5 [ T \\w\\ 2 1 + CsT sup |Mt)||g. 
Jo Jo te[o,T] 

The proof is complete thanks to the trace theorem. □ 

Lemma 2.2 (Normal trace theorem). Let w be a vector field defined on Q such that dw £ 
L 2 (fl) and divw £ L 2 (fl), and let N denote the outward unit normal vector to T. Then the 
normal trace dw ■ N exists in H (T) with the estimate 

\Bw ■ N\\ 5 < C [||flH|i a(n) + II div«;||i a(n) ] , (2.16) 

for some constant C independent of w. 

See [31] for the proof of Lemma [2~2l in the case that dw is replaced by w. For cj> £ iJ 1 (SI) , 
J T dw ■ n<j)dS = J n dw ■ D<pdx — J n div wdtfidx because we can integrate-by-parts with d 
without any boundary contributions. Thus, the identical proof given in 31j proves Lemma 
Similarly, we have 



Lemma 2.3 (Tangential trace theorem). Let w be a vector field defined on f2 such that 
dw G L 2 (fl) and curlui £ L 2 (fl), and let T\,T2 denote the unit tangent vectors to T, so that 
any vector field u onT can be uniquely written as u a T a . Then the tangential trace dw ■ r Q 
exists in H (V) with the estimate 

\dw ■ r Q | 2 _ . 5 < C \\\dw\\ 2 L2{n) + || curH|£ a(n) ] , (2.17) 

for some constant C independent of w. 

See [5] for the proof of Lemma [ 



2.4.2. An elliptic estimate which is independent of k. The following lemma is used to es- 
tablish our K-independent a priori estimates. The proof is given in Section 6 of [10] . 

Lemma 2.4. Let n > 0, s > and g £ L°°(0, T; H s (£l)) be given. Suppose that f E 
H 1 (0,T;H*(n)) satisfies 

f + Kf t =g mnx(0,T). (2.18) 

Then, 

||/IU~(o,T;J?»(a)) < Cmax{ ||/(0) || s , \\g\\L°°(o,T;H°(n))}- 

2.4.3. A technical lemma. The following technical lemma is established in |llj . 
Lemma 2.5. There exists a constant C such that 

\\dw\\n°.*<ny < C\\w\\ H o. HU) V w e H°- 5 (n). 

2.4.4. The Hodge decomposition elliptic estimates. The following Hodge-type elliptic esti- 
mate is well-known and follows from the identity —Aw = curl curl w — Ddivw, together 
with estimates divergence-form elliptic operators with Sobolev-class coefficients: 

Proposition 2.1. For an H r -class domain f2, r > 3, if w £ L 2 (il;IR 3 ) with curlw £ 
H^^-R 3 ), divw £ H 3 - 1 ^), and w ■ N\ r £ H S -^(T) for 1 < s < r, then there exists a 
constant C > depending only on Q such that 

\\w\\ s < C |H|o + || curlw|| s _i + || div io|| s _i + \dw ■ N\ s _s 

where N denotes the outward unit-normal to T. 
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In fact, the well-known version of this elliptic estimate replaces \dw-N\ s _s with |«j , iV| s _i , 
but this requires too much regularity for the unit normal N. It is easy to verify that 
having estimates for the divergence and curl of a vector field w only requires the slightly 
weaker norm \dw ■ iV| a _| estimated in order to infer full regularity of w. Also, this elliptic 
estimate is usually stated for smooth domains; the Sobolev-class regularity follows from 
the Sobolev embedding theorem and somewhat standard elliptic estimates for second-order 
elliptic operators with Sobolev-class coefficients. 

2.4.5. A polynomial-type inequality. For a constant M > 0, suppose /: f H f(t) > con- 
tinuously and satisfies 

f(t)<M + tP(f(t)), t>0, (2.19) 
where P denotes a generic polynomial function. Then for t taken sufficiently small, 

f(t) < 2M. 

3. A PARABOLIC K-APPROXIMATION OF THE SURFACE TENSION PROBLEM (TTjl 

In this section, we define a parabolic approximation of the surface tension problem (ll.7[) , 
which we term the re-problem. The K-problem is defined by adding artificial viscosity terms 
to the Euler equations and the Laplace- Young boundary condition. The salient feature of 
the K-problem is its compatibility with our energy-estimates methodology based on Proposi- 
tion [2T] in that (1) the transport structure of the Euler equations is maintained and (2) the 
momentum equations of the K-problem are equivalently expressed in the form / + «/ t = g 
with / equalling the gradient of J. These structural properties of the K-problem respectively 
yield the K-independent curl- and divergence-estimates. 

3.1. Assuming C°°-class initial data. In our construction of solutions to the surface 
tension problem (II. 7|) . we assume that the initial data (po,Uo,Q) is of C°°-class and satisfy 
the conditions ()1.13[) and (|1.14j) . as in Appendix lAl Later, in Section f6T2| we will recover 
the optimal regularity of the initial data stated in Theorem ll.il 

3.2. The parabolic approximation of the surface tension problem (|1.7p . We recall 
that an equivalent expression of the surface tension problem (|1.7p is 

vl + 2A^(p J- 1 ) lk = infix (0,T], 

p 2 J- 2 = p - ag a Pr), a0 -n onrx(0,T]. 

We have used the identity H{n) = —g a ^rj, a p -n in writing the boundary condition. 

The variables in the following problem a priori depend on the parabolic parameter k. To 
indicate this dependence, we place the symbol ~ above each of the variables. 

Definition 3.1 (The K-problem). For k > 0, we define v as the solution of 

il + 2A'y{poJ- 1 ), k ^KA'>ipoJt),k = znfix(0,T K ], (3.1a) 

p 2 J- 2 - Kpp^Jt = H K on Fx (0,T K ], (3.1b) 
(ry,u)| t=0 = (e,u ) on Q. (3.1c) 

The function H K appearing in the right-hand side of (I3.1bj) is defined as 

H K = 0(t) - o-g a ^ iaP -n - Kg af> v iafi -n, (3.2) 
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where the function (3(t) appearing in the right-hand side of (|3.2|) is defined as 

\t=o 

|t=o. (3.3) 

Remark 8. The particular artificial viscosity — KA k (poJ t ),k appearing in (|3.1a[) preserves 
the transport structure of the Euler equations. In comparison, an artificial viscosity of the 
form — kAI[A^v 1 ,k ],j would not preserve the transport structure of the Euler equations. 

Remark 9. The initial data satisfy the compatibility conditions (| 1 . 14[) so that definition p. 31) 
of j3(t) implies that (3(t) — > j3 in the limit as k tends to zero. Formally, the K-problem p. II) 
is asymptotically consistent with the surface tension problem (| 1 . T[) . 

Remark 10. For <f> G L 2 (0,T; H 1 ^)) such that • n € L 2 (0,T; H 1 ^)), the variational 
equation for the K-problem p.l[) is 

/ / Povt'4- [ [ PoJ~ 2 a%4>\ k +K [ [ PoJtafipoJ' 1 ^)^ 
Jo Jn Jo Jn Jo Jn 

+ [ [ PV)y/§<l>-n+ [ f[afj, +Kv,p] i (h i ^g a ^-n),a = O- (3.4) 
Jo Jr Jo Jr 

Since lim K ^o P(t) — f3, we have that (|3.4j) with k = is the variational equation for the 
surface tension problem ()1.7[) . Furthermore, for k > our choice of artificial viscosity 
— nA k (poJ t ),k has the nice property of not introducing any nontrivial boundary integrals in 
the variational equation p.4[) , whereas the use of the artificial viscosity — kA 3 t \A^.v 1 ,}~ for 
example, in p.4p would require additional boundary conditions to account for the tangential 
components of —kN^ AlA k v\k- 

3.3. The constant-in-time vectors v a for a = 1,2,3,4. The vector field Vt\t=o is com- 
puted using the momentum equations p.lal) . as follows: 

v t \t=o = ^KA k (poJt),k-2A k (p J^ 1 ),k^ |i=o = D(np Q divu a - 2p ). 
Similarly, for all a € N, 

d?v\ t =o = J^zr («#(po Jt), k -2A k {poJ- l ),k ) | t=0 on ft. 

This formula makes it clear that each d?v\t=o is a function of space-derivatives of the initial 
data Uq and po- We define the constant-in-time vectors v a as 

y a = ^ I ^A k (p J t ),k~2A k (p a J- 1 ), k ) \ t=Q , for a = 1,2, 3, 4. (3.5) 

Since 2A k (p J~ 1 ),k = PQ 1 a k (p l J~ 2 ),k, we have that v a — > v a , a = 1, 2, as k where v a 
are defined in Section [l~5l We use p.lbl) to compute the following identities: for a = 0, 1, 2, 3, 

d?[pp- 2 ~ *P 2 J- 3 Jt]\ t=0 = - <rg a0 V,aP -n - Kg aff v, a0 -n}\ t=0 - (3.6) 
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4. A PRIORI ESTIMATES FOR THE K-PROBLEM (I3TT 



We establish our K-independent a priori estimates in this section; the precise estimate is 
stated below in Lemma I4TT1 Our existence of solutions to the K-problem (|3.1[) is established 
in Section [5] 

For k > 0, we define the following higher-order energy function: 

E«{t) = 1 + e \\d?m\\i-a + \vm ■ mil + E ^ • "(*)ii.5-a 

rp rp 3 2 

+ I l^dVutt ■ n\l + [ ||vWi||? + E \\Kd?v(t)\\l-a + E 1^^* • »(*)l!.5-a- 

(4.1) 

Remark 11. The inequality stated in Theorem 15 . 1 1 ensures that E K (t) is continuous in time. 
We make the following definition to allow for constants to depend on 1/6: 

Definition 4.1 (Notational convention for constants depending on 1/6 > 0). We let V 

denote a generic polynomial with constant and coefficients depending on 1/8 > 0. 
We define the constant A/q > by 

ATo='P(||uo||ioo,||po||ioo). (4-2) 
We let 1Z denote generic lower- order terms satisfying 

[ 7L < N + 5 sup E K (t)+TV{ sup E K (t)). 
Jo t£[Q,T] te[o,T] 

The artificially high iJ 100 (f2)-norm in defining Ao is acceptable as the initial data (po, Uo, f2) 
is of C°°-class. We shall assume that 

^ < J < | for all t G [0, T] and at G 0. (4.3) 

The bounds (|4. 3[) are possible by taking T > sufficiently small, since thanks to Theorem l5.ll 

||J-l|U~ ( n}<C|| ! d t J\\ 2 <CVf. (4.4) 
Jo 

Lemma 4.1 (A priori estimates for the K-problem). We let v solve the n-problem (|3 - 1 [) on 
a time-interval [0,T], for some T = T K > 0. Then independent of n> 0, 



(4.5) 



sup E K {t) < [ K. 

t€[0,T] JO 

We will establish Lemma |4~T1 in the following six steps: 

Step 1: The k- independent curl-estimates. We follow [T3] in establishing 

Lemma 4.2 (The K-independent curl-estimates). 

4 ~t 3 .t 

sup V||curl9 t Q ^(t)||2_ a + / || V« curl 5«ti || o+ sup V \\ncvLT\d?v{t)\\l_ a < / K. 
* G [o,T] a=0 Jo *e[o,T] a=0 Jo 
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Proof. The Lagrangian curl of (|3.1a[) yields curl^ v t = 0. Setting 

B(A,Dv) =6.^1^,,, 
we find that (curl^ v)t = B(A, Dv). By the fundamental theorem of calculus, 

cw\r,v(t) = curlwo + / B(A,Dv). (4.6) 



Applying the gradient operator D to (14. 6|) and a second application of the fundamental 
theorem of calculus, we find that 

.Dcurl?7(t) = tDcvj\u - e.^Drf lS / A t s . 

Jo 

+ E.ji f [A t s J Df,\ s -DA s J i\ s }+ f f DB(A,Dv), (4.7) 
Jo Jo Jo 



where we have used the identity curLj Dfj = curl Dfj + e.^Drf , s J Q A t s -. 

The differentiation formula (|2.7j) equally holds for when the gradient D replaces ft; hence, 
the first three terms on the right-hand side of (|4.7[) are, with respect to the action of D 3 , 
each bounded by J Q T 1Z. 

We next analyze the highest-order term created in L f Q D 4 B(A, Dv). With (|2.7p . 

DB{A, Dv) = -e. ji [A!Dv r , q A s r v l lS +A q ]V r „ A s r Di\ s +v r , g v\ s D(Affi)], 
from which it follows that the highest-order term of D 4 B(A, Dv) is 

' [A q p%% q A s r v i , s +A q j v r , q A^D% i , s ]. 

With a relaxation of the precise structure of the summands in the integrands of the highest- 
order terms of D 4 B(A, Dv), we highlight the derivative count that results from integration 
by parts in time by writing 

pt pt pt pt pt 

I I D i B(A,Dii) = - / D 5 ij(DvAA) t + D 5 ijDvAA. 
Jo Jo Jo Jo Jo 

With such a temporal-integration-by-parts computation, the action of D 3 in (|4.7|) yields 



sup || curl^(t)|| 4 < / 1Z, 

:S[0,T1 Jo 



te[o,T] 

and by the same arguments, the action of D 3 and kD 4 in (|4.6p yield 



sup || curl€7(t)||§ + sup ||recurli;(*)||| < [ K. (4.9) 

*€[0,T] te[0,T] Jo 

By returning to the Lagrangian vorticity equation curl^ v t = 0, we find that 

cm\v t = Ej.ivi, s [ A t ], (4.10) 
Jo 

and by considering the action of D 2 and kD 3 in the identity (|4.10[) we infer that 



(4.8) 



sup || curlu f (t)||2 + sup ||recurl5 t (t)||3 < / K. (4.11) 

te[o,T] te[o,T] 
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By considering the action of Ddt, 9 t 2 , nD 2 dt, uDd 2 , and y/~Kdf in (|4.10p . and using the 
fundamental-theorem-of-calculus identity d^Vt = dfVt\t = o + J Q d^ +1 vt in the lower-order 
terms, we establish that 

SUP Y,\\ CUrW tMmi-a+Y,\\ KCUrW t dtt ^-a+ / IlKCUTlVttttWIlo < / ^. 
*e[°. r ]a=0 a=0 J JO 

(4.12) 

The sum of the inequalities (|4.8p ~ (|4.12p completes the proof. □ 

Step 2: The K-independent estimates for vutt and v t tt- We equivalently write the 
momentum equations (|3.1ap and boundary condition f|3 . lb[) of the K-problem as 

Po vl + a^(p 2 J- 2 ), k -K Po J- 1 a><(poJt),k =0 in Q x (0, T], (4.13a) 
pi J- 2 - upp-Kjt =H K on T x (0, T], (4.13b) 

where H K is given by (|3.2j) . From p. lap , we have the identity 

VfV~t = 9 7 [-2poJ -1 + np J t ]. 
Multiplying this identity by poJ^ 1 yields 

poJ~ l v t ■ fj 7 = B^[-plJ~ 2 + KplJ~ x J t } - k{po J _1 ),7 PoJt- 
We use the boundary condition (|4. 13b[) to obtain the following tangential identity for v t : 

p a j~ 1 Vffj 11 =d 1 crg^f},^ -n + Kg^v,^ -h - /3(t) - k( Po J _1 ), 7 p J t . (4-14) 
Proposition 4.1 (Energy estimates for the fourth time-differentiated problem). 

sup \\v tm (t)\\ 2 + sup ||a t 4 J(i)||o+ sup \vttfn(t)\ 2 
te[o,T] te[o,T] *e[o,T] 

+ [ T \V^dfv-n(t)\i+ l T \\^J\\l< f T K. 



Proof. Testing four time-derivatives of (|4.13al) against dfv in the L 2 (f2)-inner product, and 
integrating by parts with respect to dk in the interior integrals f Q afdf(plJ^ 2 ),k dfv 1 and 
-k J n dfdf[p J~ 1 (p Jt),k}dfv\ we find that 



dl[p vl]dfV - / dt[pU- 2 ]~atdtv\ k +n / dtlpiJ-'JtYatdtV* 
n Jn Jn 

" v ' " v ' 

+ J d?[p 2 Q J- 2 - Kf%J- 1 J t ]a$8iv i N h = K. 

v ' 

i 

It is convenient to rewrite this equation as 

\j t I P^tv? + j t I P 2 Q j- 3 \dtJ\ 2 + K f pp-i\d?J\ 2 +i = K, (4.15) 

where the identity J t = a s r v r , s implies that the error created in order to write a^dfv l ,k as 
dfj in X\ and I2 is of lower-order and so is absorbed in 1Z. 
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Rewriting the boundary integral in (|4.15|) . Using the outward normal identity (|2.9|) 
and the boundary condition (|4.13bj) with £7 = 1, we find that 

gg al3 v ttUa -dthvttt,/} -n-\ ( d t [y^g al3 } v t tt,a -nvuup -n 



99 ■<!<-■'. lit' ,j 

r z Jr 



n 



V2 

3 



+ / vtttrf -"(v^g Q ^))a vtttt ■ n - y~] ci / \HsT^ d\ff -.ap df l h j v t ttt ■ n 



4 r 

J> J d l t ~g a ^-% a0 -ft] VI 



gvtttt ■ n 



+ * jf t&«,0 {n l g aP ^) ia VffSL - ci J r d\{g af) n l )dt% a p y/§Vtm • n ■ (4.16) 



In our analysis of J T ti, i = 1, 2, 3, we adopt the convention of letting 

£ denote a function of L°°(T)-class and I a function of if°' 5 (r)-class. 
We recycle the symbols j, )a, )b, etc. in terming boundary integrals that require explanation. 

Analysis of v± in the time-integral of (|4. 1 6[) . The action of 8 a df in the tangential 
identity (|4.14[> provides that 

V t tt,a -fjry = Po 1 J[dcry (g^fj^v -fl + Kfj^V,^ -h - /9(t)) tt - lya], ( 4 -17) 

where the lower-order l ja £ H 5 (T) is given by 

2 

L fa = Kd a [(p J^ 1 ) n PoJt] tt + Vttt ■ (??, 7 Po J~ X ) iQ; + ^ C a [9°jj • (? t 3 ~ a (?7, 7 /Do J -1 )], a ■ 

0=1 

Setting fe' 3 = Pq 1 Jy/gg a ^g' yS v } s -n we use the tangential identity (|4.17[) . together with the 
outward normal differentiation formula (|2.11cl) . to find that 



ti 



= / efdyalg^fj^-nlttVtttrf-n + K / ifd ja \g» v v,^-n} tt Vttt,p-n+n. (4.18) 
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Concerning ti', we integrate by parts with respect to a time-derivative to write 

n' = - / d~iW v f)^ -n]tt d a [ifv tt , -n] T + [ I B^fj,^ -n]« 8 a [£fv tt , -n t ] 
Jr Jo Jr 



+ / / d y [r v fj^u-n}ttd a [d t efvtt >r n}+ / / B^fj,^ -n] m d a [ifv tt ,p -n] 
Jo Jr Jo Jr 



where we have further integrated by parts with respect to d a on the right-hand side. 
Interpolation and two applications of Young's inequality provide that 

|*W -n{T)\l < C\vtt ■ n{T)\ x ^Bv tt ■ n(T)|i. 5 < C\v tt • n(T)|§. 6 + 2S\Bv tt ■ h(T)\l 5 . 

Hence, 

UWlo ^ / n - 



We record that via Lemma 12.31 

\vu, a p -rir° \ 2 H -o. Hr) < C [\\B 2 v tt f + || curiaSttllg] < C\\v a \\l 

Since n t = —f) )a g' Tp v ) p -n and d 2 v t • n is in H 15 (T), we use an ff _05 (r)-duality pairing in 
the highest-order term to write 

)B< I K + C f \B 2 V t ■ n(t)|l. B |«tt,a0 -fj,* \h-° *(t) < [ K- 
Jo Jo Jo 

We use the Cauchy-Schwarz inequality for the estimate 

,T 

)c< K. 
Jo 

Regarding J Q T ) D , for the integral J Q T J r B^lg^fj,^ ■n] t tt B a £™Pv a ,p -n, we integrate by 
parts with respect to 9 7 and estimate using the Cauchy-Schwarz inequality, and for the 
integral J Q T J r d 1 \g^ 1 ' f\^ v -n\ttt £^vtt,p -n, a , we also integrate by parts with respect to <9 7 
and then estimate using an iJ -0 - 5 (r)-duality pairing. These methods work equally well 
in all terms of the spacetime-integral J J r d^ffi" f}^ v •n]ut&^Vtt,ap as well as in J Q j, 
where j is defined in (|4.18l) . Thus, we have established that 



ti < / K. (4.19) 



Analysis of L T t2 in the time-integral of (|4.16|) . We have that 



t2 = / vttt,p'nygg ap vtttt'n,a+ / vttt,p -n, a \jgg ap v U tt ■ n . (4.20) 



We write the action of df in the tangential identity (|4.14p as 

v t ttff]ii=Po 1 J B 1 (g^ v ^,^ v -n + Kg^ v v,^ v -n- p{t)) m -l 1 , (4.21a) 
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where the lower-order Z 7 € H 5 (T) is given by 

3 

f 7 = K[(p J _1 ), 7 po^t] ftt +^c a 9 t Q w • a t 4 ~'(^, 7 /JoJ _1 ) 

i=i 

Letting # = -p^ 1 Jy/gg a ^g lS rj,s a -n, we use g^HD to write 



(4.21b) 



t2a 



iZW v V,iiv ■n\tttVttt,M -n+n / d^g^vut,^ -n]vttt,i3 -n+TZ, 



where we have used integration by parts with respect to <9 7 to determine Xi a ' ■ We integrate 
by parts with respect to time in order to write 



*2a 



o Jo Jr 



o 



T 

o Jr 



ijVtt,p 7 -ng fll/ Vttt,^-n+ I TZ 



ttt - 

T 



$iv tt ,p~i -h t W v f\^ v -h]tu 



^ 7 VtUP-f Id^V^ '^IttM 



o Jr 



Letting £ = p J-\/g(p s f},>yS "ft we utilize the symmetry of £^ to exchange 9 Q and 9 7 via 
integration by parts for 

l-T 



1 



1 



K|2 



^IfiTtJa,^ -n| 



o Jr 



o Jr 



ivtt^-nvtt^uinrit-^ I I Wv tUllv -n\ z = / ft. (4.22) 



Next, integrating by parts in time, we find that 



S2cj — 



o Jr 



djig^Vu,^ -ft] n(vttt,p •«) 



K+ / j. 



o Jr 



o Jr 



W djlg^^KVtt,^ -h] \/KVtttt,0 -ft 



Interpolation and Young's inequality provide for the estimate 



Iv^dvtttt ■ n\l + CsT\dv tt ■ n\i. 5 \Kdv t t ■ n\2.5 < / ft- 
We have thus established that 



/ X2a= I 
Jo Jo 



ft. 
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The analysis of J Q T V2b is similar. We set t 1 ^ = \igg a ^g lS f),Sa -n and write 



1*2/, : / v t tt,p -n, a \Jl)g a P vttt ■ n 



r 

i l13 v t ttt,p -Vn vm ■ n = / TZ+ / x 2 b ■ 



Regarding t2b', we use the identity 

vtttt,fi -Vn = P^Adrsd^V^u -n + ng^v,^ -n - 0(t)) m - (4.23) 

where = v tt tt • fj n p poJ^ 1 + v tt tt ■ Vn (PoJ^ 1 )^ +dpL n with f 7 given by (|4.21b|) . We 
integrate by parts with respect to both dp and <9 7 in J r p^ 1 Jdp 1 \g tiV fi, llv -fi\ ttt £ l13 v t tt ■ 
where the highest-order term produced by d^-integration by parts is (|4.22[) . To estimate 
the integral where dpl-y appears, we have the choice of integration by parts with respect to 
dp or an i/~ a5 (r)-duality pairing. In the integral where vtttt • fjn {poJ~ l ),fi appears, we 
use the identity (]4.21a[) . Hence, we conclude that 



T P T r T 

t26' = 



K- / i^dpl^ut 



K 



Thus, we have established that 



r 2 < / K. (4.24) 



Analysis of f Q X3 in the time-integral of (|4. 16|) . We have that 

t3= / vut, f3 ■niVdg" ), a vtttt -n-Vc; / y/jjg aP fyf? , Q p d^ l fi 
Jr , n Jr 



^ vtttt ■ n . 



We will first establish that f$ *3a cancels with a term arising in J Q T v^b.o- Recalling that 
fit = — g lS [v 7 s ■n\v,-y, the highest-order term of hutt is — g lS [vttt,s •n]fj,-y- It follows that 

^36,0= / / g a Pg jS V,afj -V^vttt^ -nvtut ■ n+ / TZ, (4.25) 



where the estimation of the lower-order terms of requires integration by parts with 

respect to a time-derivative of v U tt- Using the differentiation formulas (|2.10a[) and (|2. 101o[) . 
we write 



) 



vttt-.fi ■n(s/§g a0 ),a vtttt ■ n . 



Thus, 

i-T r T pT 

K. 



I ^3a~ t36,0 = / 

Jo Jo Jo 
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The terms t3b,i and X3b,2 are analyzed by integrating by parts with respect to a time- 
derivative of vutt, followed by elementary estimates. Thus, 

r.36,1 + / t 3 6 >2 = / TZ. 



We next examine x^b,3 — J r Vgg a ^v t t, a p 'fit vtttt ■ After integration by parts with 
respect to a time-derivative of Vutt, we find that 



Letting = \fgg a ^g' y5 v,s -n, we have that 



r 



*zb,3 = / QpVttup-rj^Vtua-n - / / v t tt,pvl t , a -(rV ff ,p P aj3 ) t 



£lpVtttt,p ■fj n v tt , a -n = - )+ TZ 



i 

By using the identity (|4.23|) in f$ }, we note that the term corresponding with Z 7 has an 
elementary estimate after integration by parts with respect to dp . The remaining terms are 
similarly analyzed, except for J Q T J r 6gPo 1 Jdp^g^vtt^ii -n] vtt,a -n in which we integrate 
by parts with respect to both dp and <9 7 . For a certain highest-order term created by 
<9^ 7 -integration by parts, we form an exact derivative: setting £~ = Pq 1 J\/gg lS v,s -n, 

£j [g al3 v u ,ap ■n], 1 g ,ll 'vu,^u -n = -~ J J d 1 i 1 \g >lv vtu^v -n\ 2 . 
This establishes that f Q T j = J Q T 7?.. We thus conclude that 

T rT 

t3fc,3 = TZ. 

Jo Jo 
Hence, 

/ t 3 < / ft. (4.26) 
Jo Jo 

Analysis of Jg T t4 in the time-integral of (14.16|) . We integrate by parts with respect to 
a time derivative of Vut in fT V4 and, if need be, spatially integrate by parts. For example, 
letting t4 = X)f=i r 4Ji we nn d that after integration by parts with respect to time, 



(V§g al3 )t('n,ap •h)ttttvttt ■ n+ TZ 
'0 Jo Jr Jo 

r T r T 

v t tt, a -dp[n{^gg al3 )tvttt-n]+ 11= TZ. 
10 Jr Jo Jo 

Similarly, integration by parts with respect to time provides for the expression 

V4,2 + / t4,3 = / TZ. 
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Finally, using the differentiation formulas (|2.10ap and (|2.10b|) . 



T 

o Jo Jr 



^/§(2§ a ^ P -§ aP g^)vm^'n,uf}, a p-hv t m-h+ 7Z = 71 



where the second equality follows from our above analysis of J Q t2b- 
Hence, 



T 

71. (4.27) 



Jo Jo 

Analysis of Jq ts in the time-integral of (|4.16p . In the first term defining t$, we 

integrate by parts with respect to B a when dp acts on nP . This term is bounded by TZ, 
since thanks to Lemma |2. 11 



\VKVtHt\o<C / \VKVtm\o.25 < s / \\\^vtttt\\l + C S T sup \\i utt (t)\\l. 
Jo Jo te[o,T] 

The remaining terms are similarly estimated. Hence, 

t 5 < [ 71. (4.28) 



Rewriting the equation (|4.15p . Summing the inequalities ([3~T5]) . (|4~24|) . (|4T26|) . (|4T27)) . 
(j4T28| yields 

£ /%< ( 4 - 29 ) 
Thanks to the inequality (I4.29|) and the identity (|4.16l) . we equivalently write (14. 15[) as 

+ / ^Ktt.r^f^tttt.a'H / Po^lv^Jl 2 = ft. (4.30) 
The time-integral of (|4.30[) completes the proof. □ 
Via Proposition the estimates of Lemma 14.21 and Proposition 14.11 imply the following 
Proposition 4.2 (The K-independent estimates for t>ut and y/KVutt)- 



sup ||w«t(*)||i+ / WVkvuuWi < / ^ 
te[o,T] Jo Jo 

Proof. The fundamental theorem of calculus provides for 

divd?f, = a s r d?7r, s -d?fi r , s f a t 8 r , a = 0,1,2,3,4. (4.31) 

Jo 

The identity Jt = a s r v r , s implies that for a — 1,2,3, dfJt is equal to a s r d^v r lS +j, where j 
scales like Dd^~ l v. The estimate for dfj stated in Proposition 14.11 therefore implies that 

sup \\divd m (t)\\ 2 < [ 71. (4.32) 

t€[0,T] JO 



24 D. COUTAND, J. HOLE, AND S. SHKOLLER 

Using the identity N = n — J Q dth, and the estimate 

\d%v- f d t h\l b _ a <TC\\d?v\\l_ a for a = 0,1, 2, 3, 
Jo 

we infer from the trace-estimate stated in Proposition 14.11 that 

sup \Bv ttt (t) ■ N\\ 5 < [ n. 
te[o,T] Jo 

Thanks to the curl-estimates of Lemma 14.21 Proposition 12.11 therefore establishes that 

sup ||«t«(t)lli < / 7^■ 

te[o,T] Jo 

A similar analysis establishes the L 2 (0, T; H 1 (fi))-estimate for y/nvutt- D 

We now improve the normal trace-estimates of Proposition 14.11 
Proposition 4.3 (An improved K-independent normal trace-estimate for i>ttt)- 

sup \v ttt -fi{t)\l < I n. 

te[o,T] Jo 
Proof. The estimate of Proposition 14.21 together with the trace theorem implies that 

f T 

sup l^tttlo 5 < R~ 
te[o,T] Jo 

Combining this estimate with the normal trace-estimate of Proposition 14.11 completes the 
proof. □ 

Step 3: The K-independent divergence- and normal trace-estimates. We equiva- 
lcntly write the approximate surface tension problem (|4.13a[) as 

2 Po j~ 2 A k J, k +Kp A k J t , k = v\ + [2J- 1 - K Jt]A k p , k , 

or equivalently, setting k! = k/2, 

J- 2 A k J, k +K'A k J t , k = ^-[vt + pj- 1 ~Kj t ]A k Po , k }. (4.33) 

Vt 

The fundamental theorem of calculus provides that (I4.33[) is equivalently written as 

DJ + n'DJt = V t - [J,k f dt{.r 2 A k ) + nJ t ,k I d t A k ] . (4.34) 
Jo Jo 

3 

Lemma 4.3 (Estimates for J tu and K Jtttt via Lemma |2~4"|) . 

2 , ,„,„ II „ j /j-M|2 



SUp \\J ttt (t)\\{ + SUp \\KJtttt(t)\\t< / n. 

te[o,T] te[o,T] Jo 

Proof. Taking three time-derivatives of (|4.34[) produces an equation of the form / ' + nf t = g: 

DJ m + Kd t (DJ m ) = d?V t + jut- (4.35) 

By (|4.33p we have that dfVt scales like v U tt + T Dvttt + nd^J. Proposition 14.21 therefore 
implies that \\dfV t (t)\\l < J* K. According to (|4~51|) . we have that jut scales like T DJ U t + 
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T nDJutt + T Dvttt- So, ||Jttt(i)||§ < J Q T ' TZ. The fundamental theorem of calculus provides 
a good estimate for J U t in L°°(0, T; L 2 (f2)). Thus, Lemma li?T4l implies that 

T 



sup \\J t tt{t)\\l < f n 

sefo.Tl Jo 



te[o,T] 

We similarly infer from equation (|4.35l) and the fundamental theorem of calculus that 



sup \\nJ tttt {t)\\l < f K. 
;efo.Tl Jo 



te[o,T] 

This completes the proof. □ 
Lemma 4.4 (Normal trace-estimates for v t t and KVttt)- 



sup \dv tt ■ n(t)\ .5 + SU P l^ttt • ^(i)| .5 < / ^- 
t€[0,T] te[o,T] Jo 

Proof. The fundamental theorem of calculus implies the desired estimates. For example, 

\dv tt ■ n(t)\l 5 < C\v tt ■ n(t)\o, 5 \dv a ■ n{i)\i.s 

<Mo + C s f \\vttt\\i+S sup \8v tt ■ n(t)\l 5 < [ U. 
Jo te[o,T] Jo 

□ 

Step 4: The K-independent higher-order estimates via Proposition 12. li The 

divergence- and normal trace-estimates obtained in Step 3, together with the curl-estimates 
of Lemma 14. 2\ imply a good estimate for vtt ■ Hence, successively repeating Step 3 yields 

2 2 J 7 

sup J2\\d?v(t)\\La+ sup J2\\Kd?v t (t)\\i_ a < [ K. (4.36) 
te[o,T] a=0 *e[o,T] Q=0 Jo 

We recall the boundary condition (|3.1b|) is 

~g aP V, a p -h + Kg al3 v, af) -n = p(t) - p 2 J- 2 + Kpp^Jt . (4.37) 

V V ' 

J 

Writing the boundary condition (|4.37|) as an equation of the form / + nf t = g yields 

g a0 V,af) -n + Kd t (g af, V,*f) ■«) =3+ nH, a p -{ng a %. (4.38) 

As nf) is of the same regularity as ry, the flow map for the K-problem does not witness an 
improved boundary-regularity. The fundamental theorem of calculus, however, does imply 
a good estimate for the right-hand side of (|4. 38[) . Hence, we infer from Lemma HOI that the 
normal trace of fj and kv each has a good estimate in L°°(0, T; iJ 4 5 (r)): 

sup \\fj(t)\\l+ sup \\Kv(t)\\ 2 5 < / K. (4.39) 
te[o,T] *e[o,T] Jo 

We collect the estimates (|4.36|) and (|4.39|) in the following 
Proposition 4.4 (Higher-order n- independent estimates). 

3 3 rp 

sup Y i mml-a+ ^ P En^(*)ii5-a< / 

tG[0,T] Q=0 te[0,T] Q=Q Jo 
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Step 5: The K-independent improved boundary-regularity estimates. With the 
estimates provided by Proposition 14.41 for ||<9°?7||5_ a and ||K<9 t a £||5_ a , a = 0, 1, 2, 3, we are in 
position to establish the improved boundary-regularity estimates. 

Proposition 4.5 (The K-independent improved boundary-regularity via Lemma |2.4[) . 

2 2 x 

sup J2\d 2 d?v-m\i. 5 - a + sup J2\B 2 d?v t -m\i. 5 - a < I k. 

te[o,T] Q=0 te[0,r] o=0 Jo 

Proof. Taking three time-derivatives of (|4.37[) . we obtain 

g a/3 v tt ,ai3 -n + ndt{g al3 vtt,af3 ■n) = j - [df (g a/3 77 >a p ■n) - g a ^v tt , a p •»] 

v v ' 

3A 

- [Kdl(g aff v, a p -n) - Kd t (g afi v tt , a p -ft)] . (4.40) 



3B 

By Proposition H31 the right-hand side of flOD} has a good estimate in L°°(0, T\ H 5 (T)). 
Lemma 12.41 therefore provides that 

sup \5 2 v tt ■ n(t)\l 5 < f n. 
te[o,T] Jo 

Since nd t (g al3 vtt,ai3 -n) = ng al3 v t tt, a p -n + KV tt ,aP •("■3 Q/3 )t, it follows from (|4.40j) that 



sup \nd 2 v tu -«(t)lo5 < / n 

:G[0,T1 JO 



te[o,T] 

The higher-order estimates are similarly established. □ 

Step 6: Concluding the proof of Lemma 14. li The sum of the estimates given in 
Propositions I4.1H4.5I competes the proof of Lemma 14.11 Taking 5 sufficiently small in the 
inequality (|4.5|) yields a polynomial- type inequality of the form (|2.19j) . Hence, for sufficiently 
small T > and independent of k > 0, 

sup E K (t) < 2j\f , (4.41) 

te[o,T] 



where the higher-order energy function E K (t) is defined in (14.1[) . 

5. Construction of solutions to the k-problem (|3.1[) 

In this section, we prove the following 

Theorem 5.1 (Solutions to the K-problem). For C°° -class initial data (po, uo, £1) satisfying 
the conditions (|1.13jl and (| 1 . 14(1 . and for some T — T R > 0, there exists a unique solution v 
to the K-problem (|3.ip verifying (v, Vt, . . . , Vtttt)\t=o — (^o, v i, ■ • ■ , v 4) an d 

sup || W*)llo + f T \vttu -n\ 2 +J2 I* \\d?nl- 2a < oo. (5.1) 

*6[0,T] Jo „- n Jo 



a=0 



Remark 12. We recall that the initial data v a , a = 1, 2,3,4, is defined in Section 1373 



We establish Theorem 15.11 via a succession of two asymptotic estimates that correspond 
with two intermediate approximate problems defined below. Each of the two intermediate 
problems involves the use of a convolution operator that smooths in the direction tangent 
to the moving boundary. We use e > and fj, > as the smoothing parameters. 
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Our first intermediate problem, which we call the Ke-problem, is defined by smoothing 
the moving domain of the K-problem. The overall structure of the Ke-problem matches 
that of the ^-problem. It naturally follows that the e-independent a priori estimates closely 
resemble the K-independent a priori estimates of Section [4] The e = formal limit of the 
Ke-problem is the K-problem. 

Our second intermediate problem, which we call the /i-problem, is defined by smoothing 
the re-artificial viscosity term Kg"^v, a p -n e appearing in the boundary condition of the «e- 
problem. The /i-problem is a nonlinear heat-type problem with Neumann-type boundary 
conditions. The /i — formal limit of the heat-type /i-problem is equivalent to the Ke- 
problem. The key to obtaining the /i-independent a priori estimates is that the diffusive 
term of the heat-type /i-problem yields a trace-estimate on the boundary T. 

5.1. Horizontal convolution- by- layers. For e > 0, we let < <p e £ C^°(M. 2 ) with 
spt(<£> £ ) C 5(0, e) denote the family of standard mollifiers on M 2 . With xn — (xi,X2), 
we define the operation of horizontal convolution-by-layers as follows: 

Aef(x h ,x 3 ) = / (pe(x h - yh)f(yh:X 3 )dy h for / e Lj oc (R 2 ). 
Jr 2 

By standard properties of convolution, there exists a constant C which is independent of e, 
such that for s > 0, 

|AJ| S <C|/| S y.feH s (T). 

Furthermore, 

e|A £ /|! < C\f\ V/eL 2 (r). (5.2) 
We recall the local coordinates near T are defined in Section [5T3] We set 

K 

e = mindist(spt& o 9 h dB+ \ T>i). (5.3) 
We define the horizontally-convolved vector field v e on T of a given vector v by 

K 

v e = J]A e [te U )o^]o^ 1 onT. (5.4) 
i=i 

Given a sufficiently smooth vector field v, we set fj = e + J* v in il and fj e = e + f v e on T. 
We define ( e to be the solution of the following time-dependent elliptic Dirichlet problem: 

AC = At? in O, (5.5a) 

C = r?e on T. (5.5b) 

We define the following e-approximate Lagrangian variables: 

;4 e = [-DCe]~\ Je = det£>Ce, fl e = JeA, [g e ] a = Ce, a <ei/3 , and y^g\h e = [d e ] T N. 

5.2. The Ke-problem and its a priori estimates. We define an intermediate approxi- 
mate problem, which we will refer to as the Ke-problem, that is asymptotically consistent 
with the K-problem (|3.ip . To indicate the dependence on the smoothing parameter e of all 
the variables in the Ke-problem, we place the symbol " over each of the variables. In the ree- 
problem, we smooth the moving boundary and use the corresponding twice-mean-curvature 
function 
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Definition 5.1 (The Ke-problem) . For k > and e > 0, we define v as the solution of 

v\ + 2[A e ] k i {p Q J- x ), k -k[A$ (po&U = in Si x (0,T K (e)], (5.6a) 

p\ J- 2 - nplJ- x J t =H e on T x (0, T K (e)}, (5.6b) 

(Ce, ^)j*=o = (e,uo) on SI. (5.6c) 

The function St appearing in the left-hand side of equations (|5.6al) and (|5.6b[) is defined as 

Jt = j € div^v. (5.7) 
The function H e appearing in the right-hand side of (|5.6bj) is defined as 

H e = /3 e (t) - agf{ e , a0 -h e - ngfv iap ■n e . (5.8) 
The function /3 e (t) appearing in the right-hand side of (|5.8|l is defined as 

3 t a 

3 

+ K E ^[-plJ-'Jt+gfv,^ -h e ]\ t=0 . (5.9) 

Remark 13. Elliptic regularity provides that The e-approximate Lagrangian flow map £ € 
defined by (|5.5[) satisfies ||Ce||s < CII^IU f° r a positive constant C independent of e. Standard 
properties of convolution provide that setting e = in (|5.5[) yields £ e = fj. It follows that 
the Ke-problem (|5.6I) is asymptotically consistent with the K-problem (|3.1[) . 

5.2.1. T/ie constant-in-time vectors v Q /or a = 1,2,3,4. Letting w solve (I5.6p . the vector 
field d^v\ t=0 for all a g N is computed as follows: 

a?*lt=o = J^zr (^[A e ] fc (po^ t ),fe -2[A e ] fc (p Jr 1 ).* ) U on n - 

This formula makes it clear that each dfV\ t= o is a function of space-derivatives of the initial 
data uq, A e uo and po- We define the constant-in-time vectors v a as 

v - = («[i £ ] fc (poJt),fe-2[ie] fc (p Jr 1 )>fc) for a =1,2, 3, 4. (5.10) 

We have that v a — > v a as e — > where v a are defined in Section 13.31 We use (I5.6b|) and the 
definition (|5.9p of /3 e (i) to compute the following identities: for a = 0, 1, 2, 3, 

S^Jf 2 - /vgJf \£]| t=0 = d?lf3e(t)}\ t=0 - dfto-gfLap -n t + K 9 fv, a& -h t ]\ t=0 . (5.11) 

5.2.2. A priori estimates for the ne-problem (|5.6p . For e > 0, we define the following higher- 
order energy function: 

E*(t) = l + \\v tt (t)\\l+ tlVtfflelt + Yl t P?v\\l-2a- (5-12) 
J a=0 J ° 

Definition 5.2 (Notational convention for constants depending on 1/Sk > 0). We let V 
denote a generic polynomial with constant and coefficients depending on 1/5k > 0. 
We define the constant A/"o > by 

Afc =P(|KI|iooJpo||ioo)- 
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We let 1Z denote generic lower-order terms satisfying 

I iL<jif a +8 sup E £ {t)+TV( sup E e (t)). 
Jo te[o,T] te[o,T] 

We assume T > is taken sufficiently small to ensure that 

1 3 i 3 

- < J < - and - < J e < - for all f € [0, T] and x € fl 

2 ~ - 2 2 ~ 2 

Lemma 5.1 (A priori estimates for the Ke-problem). We let v solve the ne-problem ([5 
on a time-interval [0,T] for some T — T K (e) > 0. Then independent of e, 

sup E e {t) < [ Tl. (5.13) 

te[o,T] Jo 

We will establish Lemma \5. II in the following four steps: 

Step 1: The e-independent curl-estimates. Taking the e-approximate Lagrangian curl 
of the equations ()5.6a|) yields 

curl^ v t = 0. (5.14) 
Integrating the identity (15. 14|) in time from to t £ (0,T] provides that 

curl^ v = curl m + / 9 t [A e ]^\ s . (5.15) 
Jo 

We may therefore infer the curl estimates from Lemma 14.21 We record this fact as 
Lemma 5.2 (The e-independent curl-estimates for d a v, for a = 0, 1, 2). 



J2 [ T \\cUT\d?v\\l-2a< f T n. 
nJO J0 



Step 2: The e-independent estimate for vu- We equivalently write the momentum 
equations (|5.6a[) as 

povl + [a^{p 2 J- 2 ) lk -KpoJ-^a^ipoJtU = in Q X (0,T K (e)]. (5.16) 

Lemma 5.3 (Energy estimates for the action of df in (15.161) ). 



sup ||u«(t)|lo+ / \\Jttt\\o+ / \dv tt ■ n e \ z < / U. 

t€[0,T] Jo Jo Jo 

Proof. Testing the action of df in (|5 . 1 6[) against v u in the L 2 (17)-inner product and in- 
tegrating by parts with respect to dk in the interior integrals f n {d e ) k 'df(p 2 ] J~ 2 ),k v\ t and 
~ K In N ? d t [poJ^iPoJt)^] v\ t yields 



n 



J d 2 [ P lJ- 2 - np 2 J^J t ][a t ] k iV l t N k = K. (5.17) 



IT 

i 
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It is convenient to write the equation (|5.17|) as 



~ / M'-nf + K I i>U 1 j,uf + i = n. (-'.IN) 



where the error created in order to write [d c ]f dfv l ,k as Jut hi % is of lower-order and so is 
absorbed in 1Z. 

We write the boundary integral i appearing in the left-hand side of (15.18[) as 



H t tW geV tt ■ n e , (5.19) 

where we have used the boundary condition (|5.6b[) and the formula (|2.9[) . The definition 
provides that H e = (3 t (t) — erg" 13 Caap m fi>e — Kg"^v, a /3 -h e . Thus, 



« / y/lu9?Pvtt,p-neVtt, a -ne + k / v\ t ,p {n\y/g~ e g^ P n{) , a v{ t +TZ. (5.20) 



Thanks to Lemma [2711 

/ \vtt\l<C j \v tt \l 25 < I ± (5.21) 
Jo Jo Jo 

Employing the Cauchy-Schwarz inequality and Young's inequality with S > thus yields 

cT pT pT 

) - \ v tt ,p-n eia ^ e g^vtt-n e < I 1Z. 
Jo Jo 



V 

Integrating by parts with respect to dp, we similarly have that L )' < f Q 1Z. We have 
therefore established that the identity (|5.20j) is equivalently written as 

i = K J VdeSt^fitt'P ■ fl cVtt:a -fie + 1Z. 

Using this identity for i in the time-integral of (|5.18[) completes the proof. □ 
Proposition 5.1 (The e-independent estimates for v tt ). 



sup ||««(*)||o+ / \vtt-n e \i+ ||u«||i<^. 
te[o,T] Jo Jo 

Proof. The desired L°°(0, T; L 2 (f2))-estimate is provided by Lemma 15.31 The inequality 
(15.211) and the trace-estimate stated in Lemma [5.31 establish the desired L 2 (0, T; iJ 1 (r))- 
estimate. We infer from the arguments proving Proposition 14.21 that Lemma 15.31 implies a 
divergence- and normal trace-estimate for v t t- Thanks to the curl-estimate for v t t stated in 
Lemma I5~2l we have by Proposition 12.11 that the proof is complete. □ 

Step 3: The e-independent estimate for Vf 
Proposition 5.2 (The e-independent estimate for vt). 



/ INI§< / n. 

Jo Jo 
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Proof. We write the boundary condition (|5.6b|) as 

Kg? P v, a i3 -n e = KplJ~ x J t - ppr 2 + Pe{t) - ag^ ( e , a -h e . 



Taking a time-derivative of this equation yields 

Kg^vt, a /3 -n e =j t + Kv, a p •a t (n e g e Q/5 ). (5.22) 
The right-hand side of (|5.22[) scales like Td 2 v t - We infer that 

52- - 1 2 



\d z v t ■ h e \ 2 5 < / K. (5.23) 
Jo 

We equivalently write the momentum equations (I5.16|) as 

[A^J t , k =(Kp )- 1 [vi + 2[A e ]U P oJr 1 ),k-^JMiPo,k] ■ 

Using the fundamental theorem of calculus, we have that 

DJt = %+J t , k [ d t [A e }l (5.24) 
Jo 

Since the right-hand side of (|5.24p scales like v t + Dv, a time-derivative of (|5.24l) yields 

/ WDJttfi < f n. 

Jo Jo 
Lemma IQ1 provides a good estimate for J tt in L 2 (0, T; L 2 (f2)). Thus, 

\\Jtt\H < / ^- (5-25) 



Via Proposition 12.11 the curl-, normal trace- and divergence-estimates for v t , respectively 
given in Lemma 15.21 and inequalities (|5 . 23[) and (|5 .25[) , complete the proof. □ 

Step 4: Concluding the proof of Lemma 15.11 Repeating Step 3 establishes 

Proposition 5.3 (The e-independent estimate for v). 

m\i< C ii. 



The sum of the e-indcpendent estimates stated in Propositions 15 . 1[ 15 .21 and 1 5.31 completes 
the proof of Lemma 15.11 



As an intermediate step in proving Theorem 15. 11 we will establish the following 

Theorem 5.2 (Solutions to the Ke-problem). For C°° -class initial data (po, uq, O) satisfying 
the conditions (|1.13|) and (|1.14|l . and for some T = T K {e) > 0, there exists a unique solution 
v to the ne-problem (15. 6|) verifying (v,Vt,..., vtttt)\t=o = ( u o> Vi, . . . , V4) and 



4 'j 1 4 r p 

SU P ||W*)llo + E/ ll«9-2a+£/ 1^"^ • ^ll- 2a < OO 
te[0,T] a=Q Jo a=0 J0 



(5.26) 



Remark 14. We recall that the initial data v Q , a = 1,2, 3, 4, is defined in Section f5. 2. II 
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5.3. Deriving a heat-type problem with Neumann-type boundary conditions. We 

will derive a nonlinear heat-type problem with Neumann-type boundary conditions which 
is equivalent to the Ke-problem 



5.3.1. Rewriting the momentum equations (I5.6a|) via (|5.15|) . Setting 

g = np J e , (5-27) 

and using the definition (|5.7[) . the momentum equations (|5.6a|) are equivalently written as 

fit - g[A t ] k (dW it v), k = div^ v[A t ] k g, k -2[A £ ] k (p J- 1 ),* . (5.28) 

Given a sufficiently smooth vector fi, the identity — Afi ~ curl curl fi — Ddivv in the e- 
approximate Lagrangian variables is the identity 

~[A t \i[[A e ] k r v lk ], ] = curl^ curl f< fi - [l £ ] s (div^ v), s . (5.29) 

Using (|5.29[) . the equations (|5.28|) are equivalently written as 

fit - g[A e Y r ([A e } k v, k ),j = £>curl£(curbj e fi) + div^ v[A t ] k g, k -2[A e ] k (p J~ l ), h . (5.30) 

Thanks to the vorticity equation (|5. 15|1 , we further have that (|5.6ap is equivalent to 

v t - g[A e }i([A,} k v, k ), 3 = IC, (5.31) 

where the vector field K, appearing in the right-hand side of (|5.31jl is defined as 

£ = gcurU (curl Uo + e. Tl f d t [A e ]p\ s ) + div, v[A e ] k g, k -2[A e } k (p Q J' 1 )* . (5.32) 

Jo 

5.3.2. Deriving a Neumann-type boundary condition for v. We decompose any vector field 
(el 3 evaluated on V into tangential components and a normal component as 

Z = egfCe,p+ene, (5.33) 

where £ a = £ • C e , Q and £ 3 = £ • h e . In the special case of a flat boundary, we may use the 
standard orthogonal unit vectors e k to write (|5.33p as 

^ = Ce a + fe 3 . (5.34) 

For the special case where (|5.33[) takes the form (|5.34[) . we have that 

— = [(curlfi) x NY + v 3 , a e l a + [div fi - v a , a }N\ (5.35) 

We define the vector b as the sum 

b = &curi + &div, (5.36) 
where 6 cur i is a tangent vector and &di v is a normal vector, respectively defined as 



6cun = ^ [gfv^p +v^g: s gfC^ •ftj^.a , (5.37a) 

'J f 



b div = -^[gfv a ,p+v a {^gf)^WH{Q]n t . (5.37b) 
For the moving boundary T e (t) — C(r), the identity (|5.35[) is written as 

%W{AS{A e } k v\ k = ((curb fi) x h^Y+nldiv, fi + -^=b\ (5.38) 

ge ^ ^ y/ge 
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or equivalently, 

W[A e \i[A € ] k r v, k = v^J-^curl^ v) x n £ + v^Jf^div^ v)n e + b. (5.39) 
Thanks to (15.15|) . we have that (|5 . 39|) multiplied by g = npoJt is equivalent to 

-2/ 



QW[A e y r [A e }«v, k = h curl + ^l Kp 2 Q (div t v)n e + gb, (5.40) 

Po 

where the tangential vector field /i CU ri appearing in the right-hand side of (I5.40P is defined 
via the e-approximate Lagrangian vorticity equation (|5.15j) as 

ft-curi = g- v /^J e " 1 (curltt +e.ji / d t [A e ]jV l , s ) x h e . (5.41) 

Jo 

We write the boundary condition (I5.6b[) as 

npl div^ v = pIJ' 2 - f3 e (t) + cxg^Q, a p -n e + ^gfv, a p -h e . (5.42) 

Setting 



-Q/J _ ,„ V 9e a p 
9e 

pp: 2 - + ^ftccx!} -ne ft el (5-44) 



(5.43) 



r V5e 

ftdiv = 

PO 

we find by using the identity (|5.42l) for np\ div; v in the equation (|5.40[) that 

eJV 3 '[le]j[^Jr«,fc = ^curl + ^div + Q" P [v, a fi '^K + Q b - (5-45) 

We define the vector field h as the sum 

h = /icurl + /idiv + 0"^[w,a/3 -n e ]fie + g[&curl + &div], (5.46) 

with the vectors /i cur i, /idiv, ^curi and 6div respectively given in (I5.41j) . (|5.44|) . (|5.37al) and 
(|5.37bp . and the function g defined in (|5.27[) . We equivalently express the identity (|5.45l) as 

QW[A e ]i[A e ] k r v, k =h. (5.47) 

Remark 15. We record that the identities (|5.42[) and (|5.44l) establish the following identity: 

g^fg~J~ x div^ v = h div ■ n e + g,f p v, a p -fie- (5.48) 

5.3.3. The heat-type ne-problem with Neumann-type boundary conditions. 

Definition 5.3 (The heat-type Ke-problem). For n > and e > given, we define v as the 
solution of the nonlinear heat-type system 

v t - Q[A e }l([A e } k v, k ),j =K infix (0,T re (e)], (5.49a) 

gN ] [A e y r [A e ] s r v, s =h + c(t) on T x (0,T K (e)], (5.49b) 

(Ce,v)\ t =o = (e,uo) on Q. (5.49c) 

The bounded, nonnegative function g is defined in (|5.27[) . The vector fields K, and h are 
respectively defined in (|5.32[) and (|5.46|) . The vector field c{t) appearing in the right-hand 
side of (|5.49b|) is defined as 

2 .a 

C W = E -d?[QN j [Az]i[A e ]r*,s -h] U=o- (5.50) 

a=0 
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Remark 16. The vector fields v a , for a = 1, 2, 3, 4, defined in Section [5.2.H are equivalently 
defined using the (a — l)th time-derivative of the momentum equations (|5.49a[) evaluated at 
time t = 0. 

Remark 17. According to Appendix [Cl the heat-type Ke-problem (I5.49[) is equivalent to the 
ne- problem (I5.6[) . 

Proposition 5.4 (Solutions to the heat-type Ke-problcm (15. 49(1 ). For C°° -class initial data 
(po,uo,£l) satisfying the conditions (| 1 . 13[) and (|1.14[) . and for some T = T K (e) > 0, there 
exists a unique v G L 2 (0, T; H 9 (Q)) with dfv G L 2 (Q,T; H 9 ~ 2a {n)), for a = 1,2,3,4, and 
Vtttt G L°°(0, T; L 2 (f2)) £/ia£ solves the nonlinear heat-type ne-problem (|5.49[) on a time- 
interval [0,T] and verifies (v, v t , . . . , v tt tt)\t=o = ("o, vi, . . . , v 4 ). 

5.4. The ^-problem and its a priori estimates. We now define our second intermediate 
problem, which we term the /i-problcm. The /i-problem is a system of nonlinear heat- 
type equations that is asymptotically consistent with the heat- type Ke-problem (I5.49j) . To 
indicate the dependence on the smoothing paramater fi of all the variables in the following 
problem, we place the symbol ° over each of the variables. 

Given a sufficiently smooth vector field w, we set ij = e + J Q v in SI and t) c — e + f Q v e on T. 
We define £ e to be the solution of the following time-dependent elliptic Dirichlet problem: 



We recall that the convolution operator A^ for fi > is defined via Section 15.11 

Definition 5.4 (The /i-problem). Given n > and e > 0, for fi > we define i as the 
solution of the nonlinear heat-type system 



A( e = A?7 in fl, 
C e = f/ e on T. 



(5.51a) 
(5.51b) 





(5.52a) 
(5.52b) 
(5.52c) 



The bounded, nonnegative function g is defined as 



g = Kp J e . 



(5.53) 



The vector field K. appearing in the right-hand side of (I5.52a[) is defined as 



t = gcwl^ (curluo + e.ji / d t [A f ] s v\ s ) + div^ v[A e ] k g, k -2[A e ] k (p 




(5.54) 



o 



The vector field h^ appearing in the right-hand side of (|5.52b[) is given by 



h" = 




"i A/x[(\/fr5,a,8 -fie 



)oe l ]])oe^n e + g[b': uil + b^ v _ , 

(5.55) 



where q"^ appearing in the right-hand side of (|5.55p is defined as 
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and h cur i, hdw, b^ url , b^ iv appearing in the right-hand side of (15.55| are defined as 



curl' div 

/icurl = Q^fFJt 1 ( CUrl M + S.ji I d t [A e 



?#\ s ) x h e , 



Po 



(5.57a) 
(5.57b) 



t, = ^[E^V6(A^[(V6«-n e ).|9oft]) od-i+v-tgfgfLsp-^La, (5.57c) 
Je 1=1 

K _ nr- 

-E^(^y[^ ia ^ +v a (VEgf),p +t S H(Q] o 9 t ) o efhz. (5.57d) 



The vector field c^(i) appearing in the right-hand side of (|5.52b|) is defined as 



t" 



(5.58) 



a=0 



Remark 18. The fixed-point solution to the ^-problem (|5.52[) is established in Appendix FBI 

Remark 19. For any h defined on T, we have that h = Yli=i 6^ — Yli=i v / ^([v / ^^°^/]°^r 1 ) • 
For fj, > 0, we define the following higher-order energy function: 

2 r t 2 k rt __ 

11 < ) : Vl\\Z-2a,-Di- 

(5.59) 



a=0 J ° a=0i=l J ° 



Definition 5.5 (Notational convention for constants depending on 1/Sne > 0). We Ze£ V 
denote a generic polynomial with constant and coefficients depending on l/5ne > 0. 
We define the constant Nq > by 

^o='P(||«o||ioo,||po||ioo). (5.60) 
We let 1Z denote generic lower- order terms satisfying 

rT 

K<J\f + S sup E»(t) + TV( sup E»(t)). 

te[o,T] te[o,r] 

Lemma 5.4 (A priori estimates for the /i-problem). We let v solve the fi-problem (|5 . 52|) 
on a time-interval [0,T] /or some T — T K (e/i) > 0. Then independent of /j. 



sup E**(t) < / ?e. 

te[o,T] Jo 
We will establish Lemma HT4l in the following four steps: 



(5.61) 



Step 1: The //-independent estimates for v tt . Testing two time-derivatives of (|5.52a| 
against Vu in the L 2 (f2)-inner product and integrating by parts in the interior integral 
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I a df(Q[A e }i([Ae] k r v,k),j)vu yields 
1 d 



™ / \vu\ 2 + I g\[A] k i t u k \ 2 - I [W + c»{t)] tt v u 



= / df(lt-(g[A e Y r ), 3 [A e } k r v, k )v tt +n. (5.62) 
n 



■R 

The definition (|5.55|l of the vector field /i M provides that 



i = -E / A4(^«^-n £ )o^]^fA4(V6^-ne)°^]- / (fc + $& v ])«*« 



z=i 



We integrate by parts with respect to 9q, in i' to find that 

K 



i'=y] / A»[(V&v t t,p-foe)°0i]3e^ A »[(VTmt,a-foe)°9i] 

V / A^[(v tt ,i3 in e ^i), a ) o 9i]g^A^[(^iv tt ■ h e ) o 9i] 



Z=l ' 



K 



/ a m[(v / 6'"«)/3 ■fie) o 0i]d a g*fl A Al [(y / ^u tt • n £ ) o 0,] 



1=1 



K 



Y] / A |t [(- N /6«tt, j 9-* e )ofli]fl"fA M [(Wtt ■ (fleVTl),*) ° 



We employ an H °' 5 (2?;)-duality pairing in the integrals i' A , i' B and i' c . For example, 

I'a < C|^M ° ^|-0.5,X> ! |A /1 [(- V /^'W« • fi e ) o 9l]\ .5,V, 



< + C«|A m [(-\/6«m • he) ° ^]lo.5,x>,- 

Thanks to Lemma 12. 1[ |A M [(v^vtt ■ fo £ ) ° 9i]\q 5 ■ Dl < C\vtt\o\K[{V&Vtt ■ foe) o 9i]\i M < 
7Z + S\Afj, (in ■ foe) I? t> ■ Thus, we infer that 

K f 

i=5Z / A n[(VTivtt ■ foe),0°9i]g^iA ll [(v / l[iVtt ■ foe), a °6i] +% 

1=1 Jt >x 

where we have used the definitions (|5.57c[) and ()5.57ci[) in analyzing ). We infer from the 
time integral of (|5.62j) that 

sup ||u»(t)||§ + / IMi+y; [ t MVTAffoe)°9 l ]\lv l < [ it. (5.63) 
te[o,T] Jo l=1 Jo Jo 
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Step 2: The /i-independent estimates for vt- Similar to Step 1, testing the action of 
d 2 dt in the equations (|5.52ap against B 2 i>t in the L 2 (f2)-inner product yields 



sup \\3 2 v t (t)\\ 2 + ( T \\B 2 v t \\\ + 5^ /* \Kl(VZA ■ Ue) ° OlWi^ ^ (5.64) 
te[o,T] Jo l=1 Jo Jo 

The inequality (|5.64j) provides that \ v t \ \ 5 < it. We infer from Step 2 of the proof of 
Lemma IB . 2 1 that by viewing a time-derivative of the equations (|5.52aj) as an elliptic Dirichlet 
problem for it that 

' T \\it\\ 2 < f n- (5.65) 
Jo 

Step 3: The /i-independent estimates for v. Repeating Step 1, we test four tangential- 
derivatives of (|5.52aj) against B 4 v in the L 2 (17)-inner product and integrate by parts in the 
integral - J Q B 4 (g[A t }l ([A e ]*v, k ) :j )B 4 v to find that 

Id ' 



2 dt 



Jn l=1 JVi 

+ f g\{A e ] k B 4 v, k \ 2 = f \gy/FJr l (e.3i [ d t {A e p 4 V) „ xnj • B 4 v 



B 4 

n 



n 



K-(Q[A e y r ),j[Ae}>,k d 4 v+K. (5.66) 



We use an 77/~°- 5 (2?/)-duality pairing in analyzing the first term appearing in the right- 
hand side of (|5.66[) and conclude a good estimate thanks to the time-integral. Tangentially 
integrating by parts in I, we find that 

r T 

1< I K 



(5.67) 



Hence, the time-integral of (|5.66[) yields 

sup \\B 4 v(t)\\ 2 o+ f fd^Wl + Y. f \K[{^-h e )oe l ]\l Vi < fit. 

te[o,T] Jo ;=1 Jo Jo 

The inequality (|5.67[) yields 

1*12.6 < / & (5-68) 
o Jo 

Viewing the equations (|5.52a[) as an elliptic Dirichlet problem for v, we infer from elliptic 
regularity and the inequalities (|5.65l) and (15.68)) that 

T r T 
= 1 1 2 



v\\i< K. (5.69) 
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Step 4: Concluding the proof of Lemma 15.41 The sum of the inequalities (|5.63p . 
(|5.65jl , (|5.67|) and (|5.69|) completes the proof of Lemma 15.41 Taking S sufficiently small 
in the inequality (|5.61l) yields a polynomial- type inequality of the form (|2.19p . Hence, for 
sufficiently small T = T K (e) > and independently of [i > 0, 

sup J5"(t) < 2j\f , (5.70) 
te[o,T] 

where the higher-order energy function (t) is defined in (|5.59l) . 



5.5. The proof of Proposition l5.4l Proposition lB.ll establishes the existence and unique- 
ness of a solution v to the /i-problem ()5.52j) . Given the /i- independent estimate (|5.70p . 
standard compactness arguments provide for the existence of the strong convergence, as fi 
tends to zero, 

Ce -> ( e in L 2 (0,T;H 4 (rt)), it -> t in L 2 (0, T; H 2 (fl)), 

ii^i in L 2 (0, T; /f 4 (0)), A' 1 -> h in L 2 (0, T; H L5 (r)), 

v t -> u t in L 2 (0, T; /f 2 (0)), c^(t) -> c(t) in i 2 (0, T; if L5 (r)). 

The definitions f|5 . 54[) . (|5.55|) . (|5 . 58|) respectively define K,, h/*, c^{t). The limiting vectors K., 
h and c(t) are respectively defined by (jQ2l) . ([5T46]) and ([B30| . Letting G L 2 (0, T; iT^ft)), 
we have that the variational form of the /i-problem (|5.52p is 

{ v t -<t> + [ [A e } k r v, k (g{A e ]icl>), j = ( f K-<f>+ f ( [h» + cP(t)]-<l>. 
o Jn Jo Jn Jo Jn Jo Jr 

We infer from the strong convergence of the sequences (£ e , v, vt , K., h}\ c M (i)) that the limit 
(Cej v, Vt, fc, h, c(t)) satisfies 

Vfcf>+ f f [A t ] k r v, k (e[MW),j= f [ £-0 + r I [h + c(t)]-<t>. 
o Jn Jo Jn Jo Jn Jo Jr 

Thus, v is a solution of the nonlinear heat-type Ke-problem (|5.49l) on a time-interval [0, T] for 
some T = T K (e) > 0. Standard arguments provide that (e(0) = e and (v, v t , . . . , v U tt)\t=o = 
(uq, vi, . . . , V4). Furthermore, according to the inequality (|5.70p . 



27^ 2 
sup |M*)llo+E/ HWi-aa + E/ \d?v-n£_ 2a <2No. 

te[0,T] „-r,Jo „- n Jo 



By the higher-order regularity stated in Proposition IB. H we infer that 

4 -rp 4 „t 



sup \\vtm{t)\\l + Y, I \\9?v\\i-2a + J2 [ |^z)-n e | 2 „ 2Q <oo. (5.71) 

tefO.Tl „-nJ0 „-r>J0 



5.6. The proof of Theorem 15.21 By Lemma IC.ll the heat-type fte-problem (|5.49[) is 
equivalent to the ree- problem (|5.6p . Hence, Proposition 15.41 establishes the existence and 
uniqueness of a solution to the Ke-problem on a time-interval [0, T] for some T = T K (e) > 
verifying (v 7 i> t , ■ ■ ■ , v t ttt)\t=o — (uq, Vi, . . . , V4). The inequality (|5.71Jl establishes the 
inequality ()5.26p . 
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5.7. The proof of Theorem 15.11 A unique solution to the Ke-problem (|5.6|) exists by 
Theorem E! 

Taking S sufficiently small in the inequality f|5 . 13[) yields a polynomial-type inequality of 
the form (|2. 19[) . Hence, for sufficiently small T = T K > and independently of e > 0, 

sup E e (t) < 2j%, (5.72) 

te[a,T] 

where the higher-order energy function E e {t) is defined in (|5.12l) . 

For (j) £ L 2 (0,T;H 1 (CI)) such that (j> ■ h e G L 2 (0, T; iJ 1 (r)), the variational equation for 
the «;e-problem (|5.6[) is 

PoVft- [ [ ppr 2 [a e ]i^,k+K [ f poJtiaeliiPoJr 1 
ia Jn Jo Jn Jo Jn 

[ p e {t)VFe<P-ne+ I [[<TLp+Kv, ] i (ni^/Egf<P-n £ ), a = O. (5.73) 
lo Jr Jo Jr 

We infer from Section [5~5l and the pointwise convergence j3 e (t) — > /3(t) that the variational 

equation (|5.73l) converges to the variational equation (|3.4[) as e tends to zero. Hence, the 

e = limit v of the solutions v to the /«e-problem ()5.6[) solves the Ac-problem (|3.ip . By 

Section f5. 2. 1[ the solution v of the K-problem verifies (v,vt,..., Vtttt)\t=o — ( u o, vi, . . . , V4). 

According to (|5.72[) . 

, T 2 



sup \itt-n\i + ^2 f \\dMl-2a<^0. 

te[0,T] JO a=0 J ° 

It follows from the proof of Lemma 15.11 that by including two more time-derivatives in the 
definition (|5.12[) of the energy function E e (t), the solution v of the K-problem (|3.1| satisfies 

sup \\v mt {t)\\ 2 o+ f \vmt-n\ 2 + y2 [ \\ d ?v\\l-2 a < 00. 

:G[0,T1 JO n J0 



a=0 ' 



te[o,T] 

This establishes the inequality (|5.1I) . 

6. WELL-POSEDNESS OF THE SURFACE TENSION PROBLEM (|1.7[) 

In this section, we prove Theorem ll.ll via the K-independent a priori estimates of Section[4] 



6.1. Existence. We obtain a solution v to surface tension problem (|1.7[) in the limit of v 
as the parabolic parameter k tends to zero. According to Remark [9l f3(t) = f5 in the K = 
limit. Letting n = in Section IH we therefore conclude that the right-hand side of the 
inequality (|4.41|) depends only on M = P(E(0)). That is, for sufficiently small T > the 
energy function E(t) defined in (|1.9p satisfies 

sup E(t) < 2M Q . (6.1) 

te[o,T] 

The assumption (|4.3[) on J remains valid by taking T > even smaller if necessary. Hence, 

f(t)= Po J- 1 (t)> |A. 
Taking T > even smaller if necessary, we ensure that pit) = / o f]^ 1 ^) satisfies 

p(t) > X in Q(t). 

Since p(t) — p 2 (t) — ft > — /3 on T(t), the boundary condition (|l.lc[) establishes that 

aH(t) > -/3 on T(t). 
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6.2. Optimal regularity for the initial data. In order to obtain the i? 5 (f2)-regularity 
of our existence theory, we assumed that the given initial data is of C°°-class in Section l3~T1 
In fact, by virtue of the estimate (|6.1[) . it suffices for the regularity of the initial data to be 
such that E(Q) < oo. 

6.3. Uniqueness. We define 

E(„, t) = 1 + £ \\dMt)\\i-a + \dfv ■ n{t)\\ + J2 \ B2 d?V ■ n(t)\i s _ a . 

a=Q a=Q 

We suppose that and (7/2,^2, $2) are two solutions of the compressible surface 

tension problem (|1.7p with E(u, 0), for v = Vi,v 2 , bounded by some M > 0. 
Then by setting 

( = Vi-V2, w = v 1 -v 2 , g = [A] 2 - [/ 2 ] 2 , 
we have that (£, w, g) satisfies 

C = / w in x (0,T], (6.2a) 
Jo 

p w\ + [ax\ k iQ , k = [oa - oi]? ([/a] 2 )^ in x (0, T], (6.2b) 
p = -cr[5i] a,3 C, Q/ 3 -ni - cr[5i] a/3 ry 2 ,a/3 -ni 

+T[52] Q ^2, a /3 'n 2 on T x (0, T], (6.2c) 

CC,t",e)|t=o = (0,0,0) onfl. (6.2d) 
We will establish that w = 0. Setting 

EHt) = i2\\dK(t)\\i- a , (6.3) 

a=0 

we follow Section 2] with k set to zero. 

The Lagrangian curl operator curl^j applied to 

w\ + 2[A 1 ] k J llk =2[A 2 \ k l h lk 

provides the following vorticity equation for the difference w — v\ — V2- 

curl, u w t = 2e. Jl [A 1 ] s j ([A 2 ] k f2,k ■ (6.4) 

The time integral of the surface tension problem satisfied by w 2 yields 

v 2 - u = -2 / [A 2 ] k f2,k , 
Jo 

by which we infer that 

f T \\D 2 f[A2] k f2,k \\1<TP{ sup \\v 2 (t)\\l) < CTM . 
Jo Jo te[o.T] 

The curl-estimates for w therefore follow from the analysis proving Lemma 14.21 with the 
vorticity equation (|6.4[) replacing the homogeneous vorticity equation (|2.13p . 

Repeating the energy estimate for the fourth time-differentiated problem in Proposi- 
tion 14.11 the highest-order term of the interior forcing term J* J n d£{{a,2 — a,i] k ([f 2 ] 2 ),k ]w\ ttt 
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obeys 

/ [a2-ai]t([h} 2 )utt,kwi m < C [ \\Dvmt\\i\\wtm\\o < T P{ sup E^{t)). 
o Jfi Jo te[o,T] 

The boundary integral corresponding with 
-ct[3i] q,3 772,q/3 -ni + o[g 2 \ aSi r\2, a ii -n 2 = -v[gi} al3 r]2,ai3 • [ni - n 2 ] - <j[gi - gi\ at 3 r}2, a p -n 2 

of (|6.2c|) is similarly bounded. 
We notice that with 

\- 3 U,19, P [ n ,l s „, r - / r 2 ( 1, \-Z\nA1, 

In 



PoW [ai\ q p vutt p , q [ai] 8 r w r tttt , s = I p (Ji) [ai]«w P tt , 9 [ai]>[ ttt , 5 

Jo 

+ / pl(Ji)~ 3 [ai]lv2ttt p , q [ai]> t r ttt , s , 



we preserve an energy estimate for |j [ai]^?z;[ tt , s 1 1 § in Proposition ^. II 

Following the arguments in Step 3 of Section |4] provides control of the divergence and 
normal trace of the functions of E^ft). Proposition ^. II and the initial condition (|6.2d[) imply 
that 

sup E^t) < TP{ sup £ c (i)), 
te[o,T] te[o,T] 

for the higher-order energy function E^(t) defined in (|6.3p . Using the polynomial- type 
inequality (I2.19|) . we infer that w = as desired. 



7. The asymptotic limit as surface tension tends to zero 

In this section, we establish an existence theory for the zero surface tension limit of (|1.7I) 
via a priori estimates that are independent of the surface tension parameter a > 0. This 
asymptotic limit holds whenever the initial data satisfies the Taylor sign condition ([1.15)1 . 
and provides the following: 

< i/ jT \B\ -n\ 2 < -jT -±- N i4a k MJ- 2 ),k \d A V • n\ 2 . 

We recall that according to Theorem ll.il a solution to (|1.7[) satisfies 

sup W d tV(t)\\l- a + \vtu ■ n{t)\\ + Y, \d 2 d> ■ n(t)\l 6 _ a < C a , (7.1) 

^[0,T] a=Q Q 

for a finite bound C a > depending on 1/a. The Taylor-sign-condition assumption of 
Theorem 11.21 provides for <r- independent a priori estimates under the higher-order energy 
function E a (t) defined below in ([73]) . 

7.1. Assuming C^-class initial data. In our construction of solutions to the zero surface 
tension limit of (jl.7[) . we will assume that the initial data (po,uo,£l) is of C°°-class and 
satisfy the conditions (|1 .13)) . ([1.15)1 and ([1.16)1 . as in Appendix [A] Later, in Section [7.4. 2) 
we will recover the optimal regularity of the initial data stated in Theorem 11.21 
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7.2. The er-problem. For a > taken sufficiently small, solutions to the following problem 
are provided by Theorem ll.fi To indicate the dependence on the surface tension paramater 
a of all the variables in the following problem, we place the symbol ^ over each of the 
variables. 

Definition 7.1 (The a- problem). For a > 0, we define v as the solution of 

povt + $(p 2 J- 2 ), k = mfix(0, T a ], (7.2a) 

p 2 J- 2 ^ l3 a {t)-ag a ^V, a p-n on T X (0,T a ], (7.2b) 
(r), v)\t=o = (e,tto) on Q. (7.2c) 

The function j3 a (t) appearing in the right-hand side of (|7.2b[) is defined as 

f3 a (t) = (3 + Y / ^ d ?iP 2 oJ~ 2 - I 3 + °9 af> ^«p -n}\ t=0 - 
Remark 20. The initial data satisfy the compatibility conditions (|1.16[) . Thus, 

A,(t) = P + a^^dng^V^ -n]\ t=0 . (7.3) 

a=0 

The <T = formal limit of f} a {t) is /3. It follows that the a-problem (|7.2|) is asymptotically 
consistent with the zero surface tension limit of (11.71). 



Remark 21. We use (I7.2b[) to compute the following identities: for a = 0, . . . , 6, 

d?[plJ~ 2 ]\ t=a = d^ a {t)\ t= v - ad?[g a(3 V^ .fl]|t=o. (7.4) 
7.3. The a priori estimates for the er-problem. For a > 0, we define 



7 5 



E°(t) = i + e ii w)iiLj_i + E n^^wiiL-ia + ll^(*)ll? 

a=0 a=0 

+E II^W)||L _ , + ||V?a?*(*)||i + E |^+-fi • #t(t)ii_ a 



a=0 

3 



a=0 " a=0 

We will allow constants to depend on 1/6 > 0: 

Definition 7.2 (Notational convention for constants depending on 1/S > 0). We let V 

denote a generic polynomial with constant and coefficients depending on 1/S > 0. 
We define the constant Mo > by 

Mq = P(|KI|ioo,||po||ioo)- (7.6) 

We let 1Z denote generic lower- order terms satisfying 

[ 1Z<Mo+S sup E a (t)+TV{ sup E a (t)). 
We infer from the estimates (|4.4|) and ()7.1|) that for T > taken sufficiently small, 



1 w 3 

2 < J < 2 for a11 i e t°' T l and x e ^ ( 7Ja ) 
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Since the initial data satisfy the Taylor sign condition (|1.15|) . we also assume that 

0<v< — \=N j aia^(plJ- 2 ), k for all t € [0, T] and x G T. (7.7b) 
V9 

Lemma 7.1 (A priori estimates for the cr-problem). We let v solve the o~ -problem (|7.2p on 
a time-interval [0, T], for some T = T a > 0. Then independent of 1 >> a > 0, 

sup E\i) < [ ft. (7.8) 



te[o,T] 

We will establish Lemma \7. II in the following nine steps: 

Step 1: The cr-independent curl-estimates. We infer the following lemma from the 
arguments proving Lemma 14.21 

Lemma 7.2 (The cr-independent curl-estimates). 

7 1 
Sup X)H CUrl ^(*)HL-la+ SU P Ell^ CUrl ^ +a "Wll2 5-ia+ SU P ||V^CUrl3 t 5 l5(i)||? 

te[o,T] a=0 ' 2 te[o,T] a=0 ■ 2 te[o,T] 

3 , T 

+ sup ^||acurl9 t a r7(t)||2 5 _ ia + sup ||acurl9 t 3 i;(t)||2 < / jl. 
te[o,T] Q=0 2 te[o,T] Jo 

Step 2: The cr-independent estimates for dj J, 9 t 6 J and y/adfv ■ n. We recall that 

f = pJ- X 

is the Lagrangian density. Using the identity J -1 J t = div^ w, we have that 

d t / 2 = -2/ 2 div^, 

a 2 / 2 = -2/ 2 div^ v t - 2(f 2 A s r ) t v r , s . 

Letting the operator —2[fA 3 i dj]J~ 1 act in the Euler equations (|7.2a|) yields 

-2/ 2 div^ v t - 2fA\ [A\{pl J- 2 ), k } ,j = vl^f 2 ,, . 
Using the Euler equations (|7.2aj) to write v\ = —pQa^f 2 ,k, we infer that f 2 satisfies 

d 2 f 2 - 2fA\ [A^f 2 lk ] r3 = -p^alf 2 ^ A\f 2 ^ -2{f 2 A\) t v\, . (7.9) 

" V ' 

F 

Since F scales like DJ + Dv, it follows that d®F is in L 2 (fl). 
Similar to the tangential identity (|4.14|) . 

p J^ 1 Vff}, 1 = d 1 o-g^rj^u-n- f3 a (t) on T, (7-10) 

thanks to the Euler equations (|7.2a|) and the Laplace- Young boundary condition (|7.2bj) . 

Proposition 7.1 (Energy estimates for the action of df in the wave- type equation (|7.9jl ). 

sup \\d 7 t J{t)\\ 2 + sup ||9 t 6 J(i)|| 2 + sup l^dfv ■ h(t)\ 2 < f It. 
te[o,T] te[o,T] te[o,T] Jo 
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Proof. Testing six time-derivatives of (|7.9|) against Pq 2 J 3 dJ f 2 in the L 2 (f2)-inner product, 
and integrating by parts with respect to dj in the integral —2 /„ p$ 1 ajdf [af f 2 ,k ] ,j dj f 2 
yields 

\ ± J^Po 2 J 3 \dJf 2 \ 2 + 2 J p^df [af/ 2 , fc ] uU); J 2 ., -2 jf p^d? \a*f 2 , k ] alWdl! 2 

s / v. > 

X i 

- J / (Po 2 ^)^ t 7 2 l 2 +2^ Ci / L 9 t 6 r o:f 2 ,/). J p {] 2 .^ ) ;f 

L Jn ,_, Jn 



i=i 



l=1 Jn Jn 



n n 



+ 2 / p^ai(J-%d![aU 2 ,k]Jdlf 2 . 



n 

(7.11) 

We have used the Cauchy-Schwarz inequality to analyze all of the terms in the right-hand 
side of (THD) except for the highest-order terms of j n df{fA j i )d 2 {A k J 2 , k ), J p^ 2 J 3 d 7 t f 2 , 
where we have used an L 4 -L 4 -L 2 Holder inequality. 
Writing d t f 2 = — 2p\.J~ 3 J t , it follows that 

djf 2 = -2plJ- s d^-2^ Cl d l t (p 2 J- 3 )dJ- l J. (7.12) 

i=i 

The identity (I7.12j) provides that the equation f|7. 1 1[) multiplied by | is equivalent to 

jJ^p 2 J- 3 \dJj\ 2 +l~i^1Z. (7.13) 
Analysis of I in f|7. 1 3[) . We equivalently write 



1= f p,, 1 r/f n'jfi j' 2 ., •» f p^dt i mf 2 *%dH 



2 

'J 



We have that 

Similar to (j7.12j) . we have that 9 t 6 / 2 is equal to —2f 2 J~ 1 dfJ plus lower-order terms. Thus, 

i a = 2^ t J^p Q 1 f i \A k dU,k\ 2 + iz. 

For Xb,i, we integrate by parts with respect to a time-derivative of d\ f 2 ,j. For example, 
Jn Jo Jn Jo 
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Since 9 t 6 a scales like Ddfv, the fundamental theorem of calculus ensures a good estimate for 
the term evaluated at time t — T. The terms L Ib,h £ = 1, • • • ,5, are similarly analyzed: 

5 
1=0 

This establishes that 

I = 2^J Pb 1 WW* I' + ( 7 - 14 ) 

Rewriting the boundary integral in (|7.13l) . We use the trace of the action of df in the 
Euler equations (I7.2a[) to write i = f r dJf 2 dJv l alN^, or equivalently, 

i= ^dJf 2 y/§d 7 t V-h. 

Using the Laplace- Young boundary condition (|7.2bp . we find that 




Jr l=0 Jr 

v ' 

7 

-aVc* f d l t g^dj' l [rj, ap -n]^djv-n+ [ d\$ a (t) y/jjdJV ■ n . 
i=i Jr iL , 



(7.15) 



Analysis of Jq h in the time-integral of (|7.15|) . The action of d a d^ in the tangential 
identity (|7.10[) provides that 

dfv, a -77, 7 = p^Jid^d^ag^V,^ -n - d?[3 a (t)) - f 7a ], (7.16a) 

where l ja is such that l ia € H~ 5 (T), \fal ia € H°- 5 (T) and is given by 

4 

l ja = dfv ■ (t?, 7 PoJ- 1 )^ +J2 cj[flf«t • 3f-'M, 7 poJ- 1 )],,, . (7.16b) 

Setting = Pq 1 J\^gg a ^g' yS v 7 s -n we use the tangential identity (|7.16[) . together with the 
outward normal differentiation formula (|2.11cl) . to find that 

h = a [ ifd^dfig^f),^ -n) dfv, p ■h+ f lf^ a ^d°v, -h +TZ. (7.17) 



r 

n 
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We integrate by parts with respect to a time-derivative of <9f v,p -n to write 



r 



T 



T 



B 7 d![r u V,^ ■n]a 3 d a [^v,p -n] + / a 2 / d^^f,,^ -h] d a [lfdfv,p -n t 



f a 2 [ B^ig^V,^ -ft] d a [d t £fd*v,p -n] + [ a 2 [ B^[g^r h ^ -n] d a [&f%*,/i -n] 
o Jr Jo Jr 

v ' V v ' 

jib ill) 

Since adfv ■ n is in H 2 - 5 (T), we may take a sufficiently small so that 

^5" ^/rr\\2 / rZr<\ rz^^.lrrw l_»5,* 



\ad 5 t v ■ h{T)\l < ^C\V^d»v(T)\i. 5 \ad»v ■ n(T)\ 2 . 5 < / K. 

Jo 

Thus, the Cauchy-Schwarz inequality provides that 

IiaIo = - / r v B 1 [dtv^-h]a 2 ifBMv,p-n]\l+ f K = f K. (7.18) 
Jr Jo Jo 

Since adfv ■ n is in H 3 5 (T), we find by use of an H _0 - 5 (r)-duality pairing that 

Wb = [ B^adfv,^ -n] ifV^Bfv^ -h t + [ K. 
Jr Jo 

v 

n 

We employ the Cauchy-Schwarz inequality to conclude that 



(7.19) 



he < / (7-20) 
Jo 

We employ the Cauchy-Schwarz inequality or an i/~ a5 (r)-duality pairing to conclude that 

)i'd< [ K. (7.21) 
Jo 

The inequalities (f7TT8j) . (j7TT9l) . (f7720|) . (j7T2T]) establish that 

ji < / (7.22) 



Analysis of §1)2 in the time-integral of (|7.15|) . We equivalently write )2 as 

h = a [ d & t v,p -h^g^djv -n, a + a [ d & t v lP -h, a ^g a0 B 7 t v ■ h . (7.23) 



)2a )2b 

The action of df in the tangential identity (|7.10p yields 

djv ■ ry, 7 = p„ 1 J%% (ag^r h ^ -h - df(3 a (t)) - Z 7 ], (7.24a) 
where l 7 is such that ^fal 1 is in H 5 (T) and is given by 

5 

f 7 = • dt'iVnPoJ- 1 )- (7.24b) 
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Letting £@ = —p 1 J\fgg a ^ g lS ij,Sa -n and using the tangential identity (|7.24p . we have that 
ha = -a 2 [ [d?v, -rd% dflg^ij,^ ■n] +TZ. 



We integrate by parts with respect to a time-derivative of dfv,/3 in order to write 
ha = -<r 2 I [^d 5 t v, p -h] n d!n„u -i 



dMd 5 t v^ -n] n d?[g^V^ -n] + 



a 2 / [P 8 5 t v,p -n] n g^dfv,^ -n+ K 



a 2 l^d?v,p-n t ] n df[g^V,^-n] 



12a" 



We write 

r-T 







ha"= / o 2 \ P 1 dtv, M -nr v dtv, llv -h+ / a 2 / d b t v,p \rd% g^dfv,^ ■ 



Using integration by parts with respect to a time-derivative of dfv,^, we conclude that 
j = / Q T TZ. Letting £ — p^ 1 J\fg g lS fj n s -n we utilize the symmetry of £@ to exchange d a 
and <9 7 via integration by parts for 



1 



W" = x / £ \r V ^v^ -rip 



1 

2 Jo 



a 2 / £d b t v ia0 -nd 5 t v^\ng^) t -- / / £ t \g^ad?v,^ -n\ 2 = / K. (7.25) 



We have thus established that 



)2a 



K. 



The analysis of Jq hb is similar. We set t 1 ^ — \/gg al3 g' y5 tj,Sa -n and write 



hb = -v' 2 / dfv ■ 







T 












n 















Jo Jr 



n,a V$9 a P V^dtv ■ nL - / o 2 / dfv>, p [W, a V$9 a0 n% 3^ 



+ I a 2 I djv, -?j, 7 dfv-h= j TZ 



)2b 



Regarding Jq hb'> we use the identity 

d 7 t v, -77, 7 = p^J[d^df(g^V,^ -n - &(*)) - (7.26) 

where ^ = <9 t 7 # • ry, 7( g poJ^ 1 + djv ■ fj n (pqJ^ 1 )^ +dpl lt with Z 7 given by (|7.24b[) . We 
integrate by parts with respect to dp-y in a 3 J Q T J r p^ 1 jdp 1 df[g IMV t), fil/ -n] P 13 dfv ■ n, where 
the highest-order term produced by d^-integration by parts is (|7.25|) . To estimate the 
integral where dpl-y appears, we have the choice of integration by parts with respect to dp 
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or an H °- 5 (T)-duality pairing. In the integral where djv ■ ry, 7 (poJ ),p appears, we use 



the identity (|7.24j) . Thus, we conclude that ^ hb = Jq and 



A <- f 

Jo Jo 



K. 



(7.27) 



Analysis of jl h in the time-integral of (|7.15j) . We equivalently write j3 as 



6 



We infer from our analysis of (|4.25p that 



j36,0 



Hence, 



J 3 



= - ^cikbj +iz. 



The terms f Q T )3b,i for I = 1, . . . , 5, are analyzed by integrating by parts with respect to a 
time-derivative of djti and then using elementary estimates. Thus, 

h = ~hb,6 + Ti- 

Integration by parts with respect to a time-derivative of djv yields 



)36,6 



fg 7 g a Pd?v, a/ 3-n t d?v-n+ / U 



136,6' 



Letting V a Q — \/gg a Pg~ l v,s -n, we once again integrate by parts with respect to time: 



)3b,6 
T 



o 



a jf %&, p d*v\ a ij?, p t a0 ) t - J a jf P af} dfyp d?v, a -n = - jf j + J 



n. 



We conclude via the tangential identity (|7.26[) that ) = TZ. Hence, 



is < / n 



(7.28) 
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Analysis of Jq U in the time-integral of (|7.15p . We integrate by parts with respect to 
a time-derivative of djv in the j4-terms and, if need be, spatially integrate by parts. For 
example, letting J4 = X^=i H,U we find that after integration by parts with respect to time, 

U,i= Fa { (^g a P) t d 7 t (rj :a0 -n)d?v-?i+ F ft 

JO JT JO 

= - [ a [ 8fv, a -d [n (V§9 a0 )td?v ■ h] + [ ft = F ft. 
Jo Jr Jo Jo 

Similarly, integration by parts with respect to time provides for the expression 

1=2 J " J ° 

Finally, using the differentiation formulas (|2. 10a[) and (|2. 10b[) . 

' V§m ati 9^ - 9 aP r v ) 7?, a/3 -n djv ■ h + { ft = ft, 

Jo Jo 



Hj - 
o Jo Jr 



where the second equality follows from our above analysis of J Q fab- 
Hence, 



U < ft- (7-29) 
Jo 

Rewriting equation (|77i73"l) . The inequalities ([7722]) . (f7727)) . ([7728]) . ([7729]) provide that the 
boundary integral i expressed as f|7.15[) satisfies 

i= \jJ V^9 aP ^dfv ia -n V^d?v, -h + ft. (7.30) 

Using the identities (|7.14p and (|7.30p . we have that (|7.13l) is equivalently written as 

| jf pi J~ 3 \dJ J| 2 + 2± jf P» l f A \A k d$U \ 2 + \jJ v y/SlT" V^d?v, a -n yfity,, -n = ft. 

(7.31) 

The time-integral of (|7.31|) completes the proof. □ 

Step 3: The energy estimates for the action of d a d?~ 2a , a = 1, 2, 3, 4. We define the 
vector a 1 ^ as fj, a -a k and let the vector ai be defined as h ■ a k . Using these vector identities, 
we decompose the cofactor matrix a as 

a k = a k a g a ^Tf^ +a$n r on T. (7.32) 
Since a k N k = -Jgn by the formula (|2.9p . it follows that 

a A = —=N j aja e . (7.33) 
V9 

According to the identity (|7.33[) . the lower-bound (|7.7b|) is equivalently stated as 

0<^<-4(Po^ 2 U- (7.34) 

Proposition 7.2 (Energy estimates for the action of ddf in the Euler equations (|7.2a|l ). 

sup \\dd?v(t)\\ 2 + sup \^dd b t v ■ h{t)\\ + sup \ddjv ■ n(t)\ 2 < f ft. 
te[o,T] te[o,T] te[o,T] Jo 
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Proof. Testing the action of ddf in the Euler equations (|7.2aj) against ddfv in the L 2 (fl)- 
inner product, and integrating by parts in the integral J Q a k ddf(p 2 J~ 2 ), k ddfv 1 yields 

/ dd^povijdd^ + f dd^[p 2 J-% k dd^- I dd?{p 2 J~ 2 }a-dd?v l , k 
Jn Jn Jn 



■x I 



J ddf[p 2 J- 2 ]a k ddfv l N k =U. (7.35) 



We used the Cauchy-Schwarz inequality or an L 4 -L 4 -L 2 Holder inequality to analyze the 
lower-order terms in ([7.351) . 



Analysis of 1 in (|7.35l) . We have that 

d 
dt 



x =i / pp- 3 \dd?J\ 2 + I 8d?(p 2 j- 2 )5d?a h v\ k + I ddUpp- 2 )d t a k dd!v\ k +n. 



Q 



Integration by parts with respect to a time-derivative of df(p^J^ 2 ) yields 



T 



Tx = - \ 8 



T 



5d?(plJ- 2 )v\ k dfa^ - / / dd?(p%J- 2 )d t (ddfa k v\ k ). (7.36) 
J Jo Jn 



Since dfJ is in H 2 (£l) and the highest-order term of dfa k scales like Ddfv, the term T\a is 
bounded by J Q 1Z thanks to the fundamental theorem of calculus. Since the analysis of lis 
is similar, we have established that X\ = 1Z. Similarly, I 2 — Thus, 

X = n. (7.37) 
Analysis of % in (|7.35|) . From the differentiation formula (|2.6[) we infer that 

%= [ J _1 (a*af -ata k r )dd b t v r lS {plJ' 2 ) lk dd G t v l +U. 
Jn 

Since a s r ddfv r , s is equal to ddf J plus lower-order terms, we write 

x = - f j-XipP' 2 )* ddfv r , s atddfff + ii 

Jn 

J- l a k {plJ- 2 ), k ddfv r aiddfv\ s - f J- X a*(j%jh 2 ), h Bdfv r aiBdfv l N s +1z. 
n Jr 



X' k 
Integrating by parts with respect to a time-derivative of ddfv l , s yields 

11 x> = j j-^ipp-^^Bdf^aidd^^l 



[ T [ J~ l a k MJ'\kddfv r atddfv\ s - [ T [ J-'a^pp^Uddfv^ai^ddfv 1 ,^ 
lo Jn Jo Jn 
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where we have again used that afddfv l , s is equal to 88® J plus lower-order terms. Lemma [23 
provides that 



X" <C [ \\dd?v\\ . 5 \\dDd?v\\ H o, Hn y < c [ \\d?v\\l 5 < [ n. 
o Jo Jo Jo 



The decomposition (|7.32l) provides that 



-k = - / ^( P 2 J- 2 ), a Bd^v ■ rj,pg af> ddfv -ti- J ^J" 1 ^ J" 2 ) )fc ddffl -nBdfv -n. 

^- v ' ^ v ' 

ki k2 



We have used the identities J 1 a = A and rj, a -A k — 5*, where 5* is the Kronecker delta, 
in writing It! . Thanks to the Laplace- Young boundary condition (|7.2bl) , we have that 



We have used the identity (I7.3[) to deduce that B(3 a (t) scales like o in L°°(r) in the first 
term on the right-hand side and we have used an L 4 -L 4 -L 2 Holder inequality in the second 
term on the right-hand side. We also have that 



k2= \jtJ r ^9J-\-al{plJ- 2 ),k}mv ■ n\ 2 



+\ [_ Bt^gJ^aliplJ- 2 )^)^ ■ n\ 2 + f v^^Po^U BBfv ■ hB3?v ■ h t . 



r 



Hence, 



X=\jJ V^ _1 [-aft(Po«7~ 2 )>* W' ,() i'' ■ »!" - C7.3S, 
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Analysis of i in (|7.35j) . Using the Laplace- Young boundary condition (|7.2bj) . we find that 
i = \~ J v^9 Q/3 \faddt/v, a -n y/aBdfv^ -h- a f y^dfv, a p ■8(hg a/3 )Bdfv ■ ft 



-° J Vdg a0 d8fv, a -h t dd%v,p -ft-i J \\flg af> )t Vadd%v, a -ft ^ddfv^ -ft 



■r 



+ / <7dd*v, p -ny^g^Bdfv ■ h, a + a \ ddfv,p -h, a y/$g a0 dd$v ■ ft 



n 



5 

+ a [ ^^.nt^J^^ + Vq f g^ %tf B\~ l W B(a^Bdfv ■ ft) 
, Jr , i=o Jr 



n 



1=1 



J2d d l t g aP dt l lv, a p -n]B(a^Bdfv • ft) + / Bdf/3 a (t)Bdfv 



(7.39) 



We integrate by parts with respect to a time-derivative of Bdfv ■ ft to write 
ji = 



^8?v, a • \d{hg afs )add b t v • ft 



iT 
It) 



v, a p ■d(ng aP )y/addtV ■ ft* 



dfv„ 



>o Jr 



>o Jr 

Since adfv ■ ft € H 2 !i (T), it follows that 



[d(ng a P)^} t *dd?v • ft ;/3 + / / y/jj%i), a ■ d(ng a P)add?v ■ ft 
1 Jo Jr l 



>/3 • 



ll 



1Z. 



Similarly, integration by parts with respect to a time-derivative of dd^v,p -ft yields 



)2 



1Z. 



io Jr 



Since \foB\v ■ ft e H 2 (T) and y/adfv € -ff 2,5 (T), for the £ = 5 term of J4, we have that 
g af3 dtV, afl dtti d(a^ddfv ■ ft) = 

g aP V^B^, a p d t n j B^^gBdlv ■ H)f T - F f g a ^dfV , a0 d t n j B^^gBd^v ■ h t ) 

Jo Jr 

dtV >a p (B t ffg afj )t B(a^Bdfv ■ ft) - J J 8 aP 8$&, a p 8 t n J B(a^Bd^v ■ ft) 
g a/3 B?V, aP d t n j B(a^Bd?v • ft) = / R. 



10 Jr 

r T ... 
K 



Jr 
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The last equality follows from the analysis of Jq h- Since the analysis of j$ and the / = 
0, . . . , 4 terms of )i are similarly established, we infer that (|7.39l) is equally written as 

{= \jtj v ^ Bd t*>" ^ ^ Bd t^P * + & (7-40) 

The time-integral of (|?35j) . Using the identities ([7157]) . ([735)1 and (fTT^Uj) in the time- 
integral of the equation (|7.35[) completes the proof. □ 

Proposition 7.3 (Energy estimates for the action of d 2 df in the Euler equations (|7.2a[l ). 
sup \\d 2 d$v{t)\\l+ sup \^d 2 v tt ffi{t)\\+ sup \d 2 v ttt -n(t)\l < [ K. 

t€[0,T] *G[0,T] f£[0,T] Jo 

Proof. Testing the action of d 2 df in the Euler equations (|7.2al) against d 2 dfv in the L 2 (Vl)- 
inner product, and integrating by parts in the integral J Q a'jid 2 df(p 2 ) J~ 2 ),k 8 2 dfv l yields 

d 2 dt[p vl\d 2 dtV + [ d 2 dta1[p 2 J~ 2 Ud 2 dtv l - [ d 2 dt[plJ- 2 ]aWdtv\ k 
n Jn Jn 

3C I 



J^ 2 dt[p 2 J^nB 2 d^N k = n. 



The analysis of 3C and i follows from the proof of Proposition 17.21 Hence, 
' 1 f Po\d 2 dtv\ 2 + H [ ^gJ- l [-aU P lJ-%k}m tt -nf 



2dt./ n 2dt Js 



+ 2di 

We have that 



^ V / ^9 Q ' 3 \/^9 2 i;t t4 , a -n yfi&Vtutf -n = K+l. (7.41) 

-i = | / p 2 j-Wj| 2 + / 9 2 a t 4 (p 2 j- 2 )a 2 a t 4 a 2 fe ^, fc + / a 2 9 4 ( P 2 j- 2 )a t a^ 2 9 t 3 ^, fe +^. 

Since d 2 d*J is in if°' 5 (fi), we use Lemma \2. 5 1 to conclude that 

T = H. (7.42) 
Using (I7.42[) in the time-integral of (I7.41[) completes the proof. □ 

We infer the following two propositions from the proof of Proposition 17.31 
Proposition 7.4 (Energy estimates for the action of d 3 df in the Euler equations (|7.2a[l ). 



sup ||a 3 u«(t)||5+ sup |V^9 3 v 4 -ri(i)| 2 + sup |<9 3 v t • < / ft. 

te[o,T] *e[o,T] t£[o,T] Jo 

Proposition 7.5 (Energy estimates for the action of d 4 in the Euler equations (|7.2a[) ). 

sup \\8 4 v(t)\\ 2 + sup \^d i f ] -h{t)\ 2 + sup \B 4 fi ■ n(t)\ 2 < f ft. 
te[o,T] te[o,T] te[o,T] Jo 
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Step 4: The u- independent higher-order estimates via Proposition [2711 
Lemma 7.3 (The cr-independent lower-order estimates for d% J, a = 0, ... ,7). 

7 ,T 

El 

a=0 



7 T 

sup £||S?J(t)||L-4«< / & 

:e[0,T]„_ n - Jo 



Proof. The a = 6 and a = 7 cases are provided by Proposition 17.11 The higher-order 
estimates for a = 0, ... ,5, are established by interpolation and the fundamental theorem of 
calculus. For example, using that dfj is in _ff 1 (f2) 



sup ||9 t 5 J(i)||?. 5 < C S sup \\d?J(t)\\j+5 sup \\d?J(t)\\l< C ft. 
te[o,T] te[o,T] te[o,T] Jo 



□ 



Lemma 7.4 (The cr-independent normal trace-estimates for dfij, a = 0, . . . , 7). 

7 7 1 

sup £ \dtm ■ N\\_ La < ft. 

te[o.T] a=a 2 Jo 

Proof. By Lemma [2721 

sup • iv| 2 _ . 5 < c sup (iia^wiig + iidiv^wn 2 ) < / T - 

te[o,T] te[o,T] v 7 Jo 



thanks to the estimates stated in Proposition l7.2l and Lemma l7.3l Using the same argument, 
the L 2 (S!)-estimates for d 3 Vtt and d 4 v respectively given in Propositions 17.41 and 17.51 and 
the divergence-estimates given by Lemma 17.31 provide that 



3 T 

sup J2\d? a m-N\i 5 _ a < ( a. 

te[o,r] a=0 Jo 



Using the fundamental theorem of calculus, the normal trace-estimates stated in Proposi- 
tions 17. 2\ 17.31 and 17.51 complete the proof. □ 

Via Proposition 12. 1[ the estimates stated in Lemmas 17.21 17.31 and 17.41 establish 
Proposition 7.6 (The cr-independent estimates for d^f), a — 0, ... ,7). 

7 T 

sup Ell«)llL-i Q < / 

te[o,T] Q=0 2 Jo 

Via Proposition 12. 11 the estimates stated in Lemma \7 . 2 1 and Proposition 17. 21 establish 
Proposition 7.7 (The cr-independent estimate for y/adfv). 



\v^d?m\\i< [ t k 

Jo 
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Step 5: The a- independent higher-order estimates via interpolation. 
Proposition 7.8 (The cr-independent estimates for y/adfrj, a — 0, ... ,3). 

i-T 



sup £iiv?w)IIL-*«< / n. 

te[o,T] a=0 2 Jo 



Proof. We note that by use of interpolation and Young's inequality, 

sup Hv^WHL < C S sup ||r/(t)||i 5 +(5 sup \\ar,{t)\\l 5 . 
t£[0,T] te[o,T] fe[o,T] 



It follows that the estimates given in Proposition 17. 61 complete the proof. □ 
Step 6: The cr-independent higher-order estimates via the Euler equations (I7.2al) . 
Proposition 7.9 (The cr-independent estimates for 9° J, a = 0, ... ,5). 

Sup EH^ J »ll4.5-ia< 

te[o,T] a=0 - Jo 

Proof. We infer from the ath time-derivative of the Euler equations (|7.2a[) that 

sup \\d?j(t)\\ 2 45 _ ia <c sup \\d?A(t)\\i 5 _ ia + c sup \\d?at(t)\\i 5 -i« + 

te[o,T] *e[0,T] 2 fe[o,T] 2 Jo 

The highest-order term of d®A scales like AADd^fi. Hence, the estimates stated in Propo- 
sition [TjU complete the proof. □ 

Given Proposition I7.9[ we now establish 
Proposition 7.10 (The cr-independent estimates for y/crdfv and y/avttt via Proposition ^. ip . 



1 

sup £iiv^ a «(*)iiL-*.< / n. 



*6[°> T 1 a=0 

Proof. Via Proposition ^. 11 the estimates stated in Lemma 1731 and Propositions [73] and [T79] 
establish the estimate for ^favttt- For yfadfv, we have that 

|v^0?« • n(i)||, g < C s \d?v ■ n(t)\l 5 + 5\ad?v ■ n(t)\t 5 < [ K, 

Jo 

thanks to the estimate for dfv stated in Proposition [721 Hence, we infer via Proposition l2.ll 
the estimate for y/adfv from the estimates stated in Lemma 17.21 and Proposition 17.91 □ 

Step 7: The a- independent higher-order estimates for crdfrj, a = 0, ... ,4. 

Lemma 7.5 (The cr-independent normal trace-estimates for crd^fj, a = 0, ... ,4). 

3 fT 

sup 5>5 2 ^-n(*)lLi a + sup \<rdfv-n(t)\l 5 < / It. 
te[o,T] a=0 2 te[o,T] Jo 

Proof. The estimates for <9" J given in Proposition 17.91 and the fundamental theorem of 
calculus provide that the ath time-derivative of the Laplace- Young boundary condition 
(|7.2bp yields the desired estimates. □ 
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Proposition 7.11 (The cr-independent estimates for crd^jj, a = 0, ... ,4). 

sup J2\\ad?m\\l 5 -i a + sup \\adfv(ml < f K. 
te[o,T] 0=0 2 te[o,T] Jo 

Proof. Since a < i/ct, we infer from the proof of Proposition 17.91 that the estimates for 
y/ad^v, a = 3,4,5, stated in Propositions 17.71 and 17.101 establish the divergence-estimates 
for avt, crvtt and avttt- Similarly, the estimates for ^/ad^v, a = 1,2, stated in Proposition l7.8l 
establish the divergence-estimates for av and afj. Via Proposition 12. 1[ the estimates stated 
in Lemmas 17.21 and 17.51 therefore complete the proof. □ 

Step 8: The <r- independent improved boundary-regularity estimates. Considering 
the action of <9 t a , a = 4, 5, 6, in the Laplace- Young boundary condition (|7.2b|) . we notice that 
the estimates stated in Propositions 17.91 and 17.111 establish 

Proposition 7.12 (The cr-independent estimates for adfv ■ h, a = 3,4,5). 

i ,T 
sup y2\adf +a v-h(t)\l_ la + sup \adfv ■ n(t)\l 5 < / It. 
te[°> T ] Q= o 2 t £l°> T ~l J o 

Step 9: Concluding the proof of Lemma 17. li The sum of the estimates given in 
Propositions I7.1H7.121 completes the proof of Lemma 17.11 Taking 6 sufficiently small in the 
inequality (|7.8j) yields a polynomial- type inequality of the form (|2.19[) . Hence, there exists 
T > that is independent of a > and verifies 

sup E a {t) < 2M . (7.43) 

te[o,T] 

7.4. The proof of Theorem [TT2l 

7.4.1. Existence. Assuming that the initial data (po,uo, f2) is of C°°-class, as in Section [7T| 
we have that Mq = P(E(0)) < oo where the higher-order energy function E(t) is de- 
fined by (|1.9p . Since the difference of the compatibility conditions (|1.14[) and (|7.4p is in 
C°°(r), we repeat the proof of Theorem 11.11 to obtain the existence of a solution v to the 
a- problem (|7.2j) . The cr-independent estimate (|7.43j) and standard compactness arguments 
establish the strong convergence, as cr tends to zero, 

j) -> rj in L 2 (0,T;H 3 - 5 (fl)), 

vt^v t in L 2 (0,T;H 2 - 5 (fl)). 

Letting cj> € L 2 (0, T; iJ 1 (J7)), we have that the variational form of (|7.2[) is 

pov t '<f>- / / plJ^ 2 a k l <j)\ k + I I f3 a (t)yf$(j)-h = cr / / y/Sg^tj^u-n^-h. 
la Jn Jo Jn Jo Jr Jo Jr 

The strong convergence of the sequences (rj, Vt) and the pointwise convergence /3 a (t) — > f3 
provide that the limit (rj, «t,/3) satisfies 

/ PoVf(f>- f f p 2 Q J~ 2 a><<f>\ k + [ [ /VS0'" = O. 
o Jn Jo Jn Jo Jr 



Thus, v is a solution of the zero surface tension limit of (|1.7j) on a nonempty time-interval 
[0, T]. Standard arguments provide that v(0) = uo and 77(0) = e. Furthermore, letting a = 
in the a priori estimates in Section 1731 we conclude that the right-hand side of inequality 
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(|7.43p depends only on Mo = P(£(0)). Hence, for sufficiently small T > the higher-order 
energy function E(t) defined in ([1.1 1|) satisfies 

sup £(t) < 2M . (7.44) 

te[o,T] 

The bounds ([7.7p remain valid by taking T > even smaller if necessary. By the arguments 
in Section [6.11 the boundedness of J in assumption (|7.7a[) implies that 

p(t) > A in U(t). 

Similarly, the lower-bound (|7.7b[) provides that p(t) = f a 77^ 1 (t) satisfies 

°<"<-S onr(t). 

on(t) 

7.4.2. Optimal regularity for the initial data. In order to obtain the -ff 4 ' 5 (f2)-regularity of 
our existence theory, we assumed that the given initial data is of C°°-class in Section 17.11 
By virtue of the estimate ([7.441) , it in fact suffices for the regularity of the initial data to be 
such that £(0) < 00. 

7.4.3. Uniqueness. We infer the uniqueness of the solution to the zero surface tension limit 
of (jl.7p by repeating the arguments given in Section [6.31 
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Appendix A. Constructing C°°-class initial data 

We demonstrate how given initial data (po,uo,^t) of optimal regularity satisfying the 
conditions (|1 . 13|) and (j!.14[) (or the conditions (|1.13l) . ([1.15[) and (j!.16[l ) we are able to 
produce asymptotically consistent C°°-class data (p , vq, ft) which satisfy similar statements 
of the conditions (fl~T3[) and (fl~T4| (or the conditions <HTT3j) , (fl~15j) and ([TT6[) ). 

A.l. The C°°-class data for the surface tension problem (jl.7p . We suppose the initial 
data (poi uq, f2) is of the optimal regularity stated in Theorem 11.11 and satisfy the conditions 
(fi~13[) and (pTT4|) . Letting 

J a = d«{J- 2 )\ t=0 and H a = d?Hir))\ t = 0) (A.l) 

we will construct C°°-class data {p ,Vo,Q) which satisfy 

aH > 4A 2 - /3 for a, /3 > 0, (A.2a) 
p > 2A > in IT, (A.2b) 

p 2 J a = d?/3 + aH a onT for a = 0,1, 2, 3. (A.2c) 
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A. 1.1. Defining the C^-class initial domain For e > 0, we let < (p e £ C^°(IR 3 ) denote 
the standard family of mollifiers with spt((p e ) c B(0,e). Recalling the local coordinates 
near T defined in Section ^. 1.51 we let C R 3 denote the C°°-class domain defined by the 
C°°-class charts 

9i = (p e *[£ Bf {9i)} in Bi for 1 < I < L, 

where for Xcl 3 the operator £x denotes a Sobolev extension operator mapping H S (X) to 
H S (M. 3 ) and where spt(£/ o 0{) C Bf C Bj. We assume that e > is taken sufficiently small 
so that the collection of open set {U{\f =1 introduced in Section T2.1.5I is an open covering 
of fi. Setting T — d£l defines a C°°-class surface. We let r] denote the C°°(r)-vector 
describing the geometry of T. In other words, drj spans the tangent space of T and, by 
letting N denote the outward-pointing unit normal vector to the surface T, we have that 
N = rj,i x 77,2 /\r),i xt7, 2 |. We let g denote the surface metric induced by 77 and let H 
denote twice the mean curvature of T. Thus, 

Hq = —g a ^ri, a p -N in C°°(r). 

We assume that the given iJ 5 -class domain C K. 3 is such that (| 1 . 14[) is satisfied. Thus, 
for e > taken sufficiently small, standard properties of convolution provide that 

P + aH > 2A 2 for cr, p > 0. (A.3) 

Hence, (|A.2a[) is satisfied. 

A. 1.2. Defining the C°° -class data (g ,Vo,fl) to satsify (|A.2[) for a — 0. Given uq £ 
H A (fl), we define uo in the C°°-class domain Q via the equation 

L 

uo = Yl& u o i\ in H 4 (n). (A.4) 

1=1 

Using the operator A e defined in Section l5TTl for e > we define the vector field Vq E C°°(S7) 
as the solution of the following elliptic Dirichlet problem: 

V a - eAy = <P £ * £n(uo) in ft, (A.5a) 

K 

V = Y, A e [fjuo ° Bi] o Gl 1 on T. (A.5b) 
Given po £ H 4 (tt), we define go in the C°°-class domain fl via the equation 

L 

Qo = ^[6poo0 / ]o0- 1 in H\fl). (A.6) 
1=1 

Assuming that po > 2 A > in Q, we let A > be such that 2A 2 > 9A 2 . It follows for Qq 
defined by dA~6l) that 

g > 3A in H. (A.7) 
For fi > and e > 0, we define g Q £ C°°(J7) to be the solution of 



g - £tAp = ipP * £ n (go) in £1, 
g Q = \/~j3 + 0-H0 on T. 



(A.8a) 
(A.8b) 
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The boundary condition (|A.8b|) implies for a = that (|A.2cj) is satisfied, since (|A.8b|) is 
equivalently stated as 

qI J = p + a H on T. (A.9) 
Elliptic regularity provides for s > 2 and a positive constant C independent of p that 

IIVMCbllff«(n) ^ c (\\^* £ ^o)\\%.- H n) + \\e f HS - iir) )- 

The boundary forcing function used in defining g Q is in C°°(T). Since the function go 
appearing in the right-hand side of (|A.8a[> is in H 4 (fl), there exists a constant C which is 
independent of p verifying 

WV^QoWh^u) < c. 

The equation (IA.8aP provides that g — pAg + tpf* * £n(f?o) m O- Hence, by standard 
properties of convolution and taking p sufficiently small, we conclude that 

II 0o _ eolU^cn) < \/Mll\/M A 0ollH 2 (n) + II V M * £ u(qo) - £»o|U°°(n) < A. 
Recalling the lower-bound (|A.7I) . we conclude that g Q > 2 A > in fi. Letting e be sufficiently 
small, the boundary condition (|A.8b|) provides that 

g > 2A > in H. (A.10) 

By (lA.lOj) and (|A.3|) . the boundary condition (|A.8bj) provides that 

P + aH Q >A\ 2 forcr,/3>0. (A.ll) 

Hence, the C°°-class data (g , Vq, fi) satisfy the conditions (|A.2|) for a = 0, as verified by 
(lA~9l) . (OTTOll and fOTjl . 

A.f .3. Formal definitions. We formally set p equal to g defined by (|A.8|) . We define 

v a+1 = -2d?(A k (p J- 1 ))\ t=0 for a = 0, . . . , 5, (A.12) 
and assume that (|A.12j) yields ■Wo = Vo defined by (|A.5|) . We make the following definitions: 
j Q = div-i; a and k a = d?(A s r v r , s )| t=0 for a = 0, . . . , 6. (A.13) 
Using the notation 

[A a } s r =d?A s r \ t =a for a = 0,..., 5, (A.14) 
it follows that k a defined in (|A.13j) satisfies 

feo=Jo, k 1 = [A 1 ] s r v r , g +j 1 and fc 2 = [A 2 ] s r v r , s +2[A 1 ] s r v[, s +j 2 . (A.15) 
We also have that J a defined by (|A.1[) is equivalently given by 

J = 1, Ji = -2fe , J 2 = 4fcQ-2fei and J 3 = -8fcg + 12fc fci - 2k 2 . (A.16) 
Using (|A.15|) and (|A.16|) . the condition (|A.2c[) that p\J a = o-H a is equivalently stated as 

divv a - 1 =j a _ 1 = --[aH a -p%J-2k a -i]+d a -i, for a = 1,2, 3, (A.17) 

/o-l 

where e£ a _i represents the lower-order terms defining fc a -i- We notice that u a , a = 1,2, 
defined by (|A.12|) is equivalently given by v\ — —2Dp and v 2 — —2{Ai} h p Q ,k —2D(p j ). 
Thus, 

divi) = fo, div«i = f 1 = — 2Ap , and divt> 2 = f 2 ~ — 2Adivu - (A. 18) 
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A. 1.4. Defining vq so that the C°° -class data (g , vq, CI) satisfy (|A.2j) for a — 0, 1. Extend- 
ing TV and t q = ij, a /\rj, a |, for a — 1,2, into fi, we decompose any vector ^ £ I 3 into 
normal and tangential components as 

e = rT«+e 3 7v inn. 

It follows for sufficiently regular £ that div£ = is equivalently written as 

divC = r,a+f , 3 +£ Q divT Q +£ 3 div AT on T. (A. 19) 

S v ' 

io 

C3 



This provides the following identity for (iV • -D)£ = £ ,3 
SAT 

We define G C°°(J7) by 



o^=io- [r,a+€ a divT a +€ 3 divJV] onT. (A.20) 



v% + eA 2 v% = p 6 *£ n (ug) in ft, (A.21a) 

q£ = Vo on T, (A.21b) 

t? 3 , = on T, (A.21c) 

where the boundary forcing function ip Q is defined as 

<Po = fo- +^0 divr Q + Vl div AT], (A.22) 

with the function / appearing in the right-hand side of (1A.22[) defined by (|A.17[) . With 
Vq denoting the tangential component of the vector defined by (|A.5|) , we define the vector 
v G C°°(n) as 

v = y„T„ + vf,N in ft. (A.23) 

The boundary condition (|A.21c[) and the definition (|A.23[) imply that 

v Q = Vo on T. (A.24) 

Thanks to (|A.19|) . the boundary condition (|A.21b|) and the definition (|A.22I) yield 

div-u = / . (A.25) 

Thanks to the definition (|A.17[) of / , we equivalently write (IA.25|) as 

qI Ji = <tH 1 on T. (A.26) 

Hence, the C°°-class data (g Q , Vq, ft) satisfy the conditions (|A.2[) for a = 0, 1, as verified by 
COT]) . (TOl) . (LOOT) and (IA361) . 

A. 1.5. Defining the C x -class data (p , i>o, ft) to satisfy (| A.2|) /or a = 0, 1, 2. According to 
(IA.18P , the condition (|A.2c[) for a = 2 may be imposed by prescribing a Dirichlet boundary 
condition for Ap . We define p G C°°(f2) to be the solution of the polyharmonic problem 

p - /tiA 3 Po = * £n(oo) in ft, (A.27a) 
-2Ap = f 1 on T, (A.27b) 

~ 2 M = V " ° n r ' (A - 27C) 

p = g on T. (A.27d) 
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The boundary-forcing function f x appearing in the right-hand side of (|A.27bj) is defined as 

fi = -^+ 2k l-\.M>U- (A.28) 

Using the boundary conditions of the poly harmonic problem (|A.21[) . we may express the 
right-hand side of (IA.28[) in terms of Vo, V\ and cp , where 

V 1 = -2Dq on T. (A.29) 

The boundary-forcing function V\ appearing in (IA.27c[) is the normal component of V\ 
and the function g is the solution of the elliptic Dirichlet problem (IA.8I) . The boundary 
conditions (|A.27c|) and (|A.27d[) imply that v\ = -2Dp Q defined by (|A.12I) satisfies 

v 1 = Vi on T. (A.30) 

Polyharmonic regularity provides for s > a positive constant C independent of p that 



2 



MPollfr*+6(n) 

dp 2 



< C(||^ * £n(go)\\h { n) + II Ap llW=(r) + |^L. +4 .. (r) + IIPollff-+«(r) 

Recalling Section fA.1.2| the arguments establishing the lower-bound (jA.lOj) provide that 

p > 2A > in H. (A.31) 

Similarly, we infer from the arguments given in Section |A. 1.41 that the boundary condition 
(|A.27b|) and the identity (|A.30I) establish that 

p\Ji = aH 2 on T. (A.32) 

Hence, the C°°-class data (p ,v ,Q), where p is the solution of (|A.27[) and v is given 
by (TQ3)) . satisfy the conditions ([A~2"l) for a = 0, 1, 2, as verified by (|A~TT|) . (jA~9)l . (|A~26|) . 
(|A~3T|) and (|A~32)) . 

Remark 22. According to (|A.f 8[) . the condition (|A.2c[) for a = 3 may be imposed by prescrib- 
ing a Dirichlet boundary condition for -^Avq in a polyharmonic problem for Vq satisfying 

A. 2. The C°°-class data for the zero surface tension limit of (jl.7[) . We suppose the 
initial data (po,Uq,Q) is of the optimal regularity stated in Theorem 11.21 and satisfy the 
conditions (|1.13[) . (I1.15[) and (|1.16l) . Then setting a = in the construction of the C°°-class 
data of Section [ATT1 J a = d^(J~ 2 )\ t= o satisfies 

J = 0, J i = 0, J 2 = and J 3 = on T. 

For p, > 0, and with go defined by (|A.6|) . g solving (|A.8j) . we define p Q e C°°(fl) to be the 
solution of the polyharmonic problem 

(A.33a) 
(A.33b) 

(A.33c) 
(A.33d) 

(A.33e) 

(A.33f) 





= ip* 1 * En(go) 


in O, 


A 2 Po 


= / 3 


on r, 


dA Po 


dAg 


on r, 


dN 


dN 


-2Ap 


= fi 


on r, 


1 dp0 


= vf 


on r, 


ON 




Po 




on r. 
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where the boundary- forcing function f 3 is defined below in (|A.36j) . By Section lA. 1.51 

Pa > A > in H. (A.34) 
According to the definition (|A.12j) of v a , 

v 3 = -2[A 2 } k p , k ~4[A 1 } k (p K ) 1 k~2D(p K 1 ) in fl, 
where K a = ^(j-^J" 1 J t )) | t=0 = df( J-\A s r v r , s ))|t=o- We compute 

J 4 = 16fco ~ 48fcgfci + 12kl + 16k k 2 - 2fc 3 . 
Setting J4 = yields an equation for fc 3 . Since fe 3 = [A 2 ]^,Vq, s +3[A 2 ]f,vl, s +3[Ai]f,v 2 \ /S +j 3 , 



1 r 



3 3 



16feo - 48fcgfci + I2kl + I6k k 2 



[A 2 ] s r v r o,s +3[A 2 ] s r v[, s +3[A 1 ]° r v r 2,s 



On the other hand, the divergence of V3 = — 2{A 2 ] k p Q ,k — 4[Ai] fc [p Ko) — 2D(p Q K 1) 
yields 

divv 3 = div [ - 2[A 2 } k p , k -4[-A 1 ] fe (p 0J ftr ),fe -2Dp Kx] - 2p on K Ui -2p AK x . 
We solve for AKi to write 



AK 1 = 

2p 



- i 3 + div [ - 2[A 2 ] fc p , fc -4[A 1 ] fe (p K ) ,t -2Dp^K^\ - 2p on K Ui 



j 3 



Using that K\ = — fcp + fci and fci = [Ai]^q, s we find that 

Aj 1 = 3 3 + A*2-A([A 1 ]j:«5, a ). 

Since Ui = —2Dp implies that jj = — 2Ap we equivalently write (|A.35I) as 

A Vo = - J [3 3 + Afco - A([Ai]X.- )] on r 



(A.35) 



(A.36) 



/ 3 



The quantity may be expressed using the boundary conditions (IA.33c[) - (lA.33fj) and 

D a VQ, a = 0, 1,2,3 may be expressed using the boundary conditions of the polyharmonic 
problems defining vq. It follows from (|A.12|) and (|A.36j) that J a — dt(J~ 2 )\t=o satisfies 



J a =d?(3 on T for a = 0,1, 2, 3, 4. 



(A.37) 



r) 2 r-1—1 

With < 2v < —^77, we may take e sufficiently small so that q defined in (|A.8[) satisfies 

0<,<-go,r. 

We then have that (|A~29)) defining the vector field V x € C°°(r) is 



M 



AT on T. 



The boundary condition (|A.33cj) of the polyharmonic problem (|A.33j) defining p provides 
that 
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The conditions (|A.33e|) . (|A.34j) and (|A.37j) establish that the C°°-class initial data (p , vo, SI) 
satisfy 

Po > A > in SI, (A.38a) 
8n 2 

-^>^>0 onT, (A.38b) 

pp a =d?/3 onT for a = 0,1, 2, 3, 4. (A.38c) 
The conditions (|A.38[) are analogous statements of the conditions (|1.13|) , (|1 . 15|) and (|1.16p . 
Remark 23. According to (|A.12|) . 

div-u 4 = div [ - 2[A 3 ] k p Q , k -6[A 2 ] k (p Q K ), k ~2[A 1 } k (p K 1 ), k -2Dp K 2 ] - 2p 0)i K 2n 

-2p AK 2 . 

Since AK 2 is approximately A 2 div«o, the boundary condition J 5 = may be imposed by 
prescribing a Dirichlet boundary condition for -^A 2 Vq in a polyharmonic problem for Vq 
satisfying Vq + eA 6 Vq = Lp e * £n(ug). 

The boundary condition Jq = is obtained by defining a Dirichlet boundary condition 
for A 3 p in a polyharmonic problem for p Q satisfying p — pA 7 p = ip* 1 * £n(f?o)- 

Appendix B. Solutions to the /^-problem (|5.52[) 

In this appendix, we construct a fixed-point solution v to the /z-problem (|5 -52[) : 

Proposition B.l (Solutions to the /i-problem). For C°° -class initial data (po,uo,Sl) satis- 
fying the conditions (| 1 . 13[) and (jl - 14[) . and for some T — T K (e/i) > 0, there exists a unique 
v G L 2 (0, T; H 9 (Sl)) solving the p-problem (15.52[) on a time-interval [0,T], with d^v G 
L 2 (0,T-H 9 - 2a {Sl)) for a = 1,2,3,4, butt G L°°(0, T; L 2 (fl)) and (v, v t , . . . , £«tt)|t=o = 
(u , vi, . . . , v 4 ). 

Remark 24. We recall that the initial data v a , a = 1,2, 3, 4, is defined in Section f5. 2 .11 

B.l. The functional framework for the fixed-point scheme. To establish the exis- 
tence and uniqueness of solutions to the /z-problem (|5.52l) . we define for T > the Hilbert 
space 

X T = {v G L 2 (0,T;H 5 (SI)) | d£v G L 2 (0, T; H 5 ~ 2a (Sl)) for a = 1,2}, 
endowed with the natural Hilbert norm 

2 

IMIx T = XI W^t v \\'L^{0,T;H'-'- 2a {^))- 
a=0 

For A4 > (where the particular value of A4 is specified later), we define the closed, 
bounded, convex subset Ct{M.) of 

v G C T (M) C X T (B.l) 
to be all v G Xt that satisfy each of the following conditions: 

(XI) (v,v t ,v tt )\t=o = (uo,vi,v 2 ), and 
(X2) \\v\\ 2 Xt <M. 

Lemma B.l (Solutions to the /i-problem in Xt). For C 00 -class initial data (po,uo,Sl) 
satisfying the conditions (jl.131) and (|1 . 14|) . and for some T — T K (ep) > 0, there exists 
a unique i) G Xy that solves the ^.-problem (|5.52[) on a time-interval [0, T] and verifies 
(v,Vt,V tt )\t=0 = (wo,Vi,v 2 ). 
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B.2. Linearizing the //-problem (|5.52[) . Letting v £ Ct(-M), we set 7] = e + J Q v and 

fj e — e + J Q v e on T. We define ( e to be the solution of the following time-dependent elliptic 
Dirichlct problem: 

= At; in fl, (B.2a) 

Ce = % on T. (B.2b) 
We define the following e-approximate Lagrangian variables: 

A £ = [DQ- 1 , J e = detL>C e , a e = J e A e , [g e ] a p = ( e , a , and y/g\n e = [a e ] T N. 

We assume that T > is given such that independently of the choice of v £ Ct(-M), the 
e-approximate Lagrangian map £ e is injective for all t £ [0, T], and that 

1-31-3 — 
2 < J (t) < 2 and 2 - - 2 f ° r a11 1 G '°' T ^ and X G n ' < ~ B ' 3 ' ) 
This is possible by the inequality (|4.4p as -D € Ct(A^) satisfies ||w||x T < -A4 . 

Definition B.l (The system of linear heat-equations for v). For v £ Ct(-M), k > 0, e > 
and [i > given, we define v to be the solution of the system of linear equations 

vt - g[A e }i([A e ^v,k),j =1C infix (0,T re (e/i)], 

gNi[A e ]l[A e } s r v, s = h* + c"(i) on T x (0, T K (efj,)], 

v\t=o = Uo on Cl. 

The bounded, nonnegative function g is defined as 

g = Kp J t . 

The vector field K, appearing in the right-hand side of (|B.4aj) is given by 



(B.4a) 
(B.4b) 
(B.4c) 



K = gcurl c - £ (curlu + £-ji / d t [A e } s jV \ s ) + div fe v[A e ] k g lk -2[A e ] k (p J~ 

Jo 

The vector field h 11 appearing in the right-hand side of (|B.4b[) is given by 

K 

hf* = h cur i + hdiv + ^2 ( A M &e,l A A« [(V&V,al3 'Tie) ° 0{\ 



(B.5) 
(B.6) 



i=i 



9 ; n e + g 



curl 1 div 



where 



Po 



and the vectors h CU ii> hdi v , 6(f url and b^ iv are defined as 

h cm \ = gVgl{Je)~ 1 (curlu a + e. jt \ d t [A e ]^v\ s ) x n e , 

Jo 

h div = ^SL \ P 2 (Je)- 2 ~ Pe(t) + O-gf&M 
Po L 



"curl — J 



i/3 -Tie 



K 



5/3 -n e 



1=1 



K 



= - E V^(Ap [ VPj 1 [5e Q ^ Q ,/3 +V a (VFegfU +** H ((e)] 



o 



J l '"e< 



(B.7) 
(B.8) 

(B.9a) 
(B.9b) 

(B.9c) 
(B.9d) 
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The vector field c M (t) appearing in the right-hand side of (|B.4bj) is given by (|5.58|) . 

B.3. Implementation of the fixed-point scheme to solve the /i-problem (|5.52|) . We 

will allow constants to depend on 1/Snefj, in our fixed-point scheme. 

Definition B.2 (Notational convention for constants depending on l/Snefj, > 0). Given 
v G Ct(-M), we let V denote a generic polynomial with constant and coefficients depending 
on 1/<5k£/i > 0. We define the constant Mq > by 

JV = P(\\u \\ 1Q0 , UpoIIioo). 

We let it denote generic lower-order terms satisfying 



f n<Mo + S\\v\\x T +TP(\\v\\ 2 XT ). 
Jo 



Lemma B.2. Given v G Ct{M), k > 0, e > and // > 0, there exists a unique v G 
solving (|B.4p and satisfying \(Xl)\ Furthermore, 

IMIx T < [ T K- (B.10) 



Proof. Standard parabolic theory provides for the existence and uniqueness of v G Xj- 
solving (|B.4[) and satisfying ! (XI )| For the purpose of establishing the estimate (|B.10[) . it is 
useful to note the scaling relations 

lC~Dv+ I D 2 v and ~ V + / Dv + K^v + h e , 
Jo Jo 

where h e G C°°(T) with \h e \ s < C e |ti| s thanks to the inequality ()5.2j) . We will establish the 
inequality (jB.lOj) in the following three steps: 

Step 1: Parabolic estimates for v u - Testing two time-derivatives of the equations (|B.4aj) 
against v tt in the L 2 (f2)-inner product and integrating by parts with respect to dj in the 
integral - J Q (g[A e ]l([A e }^v , k ),j )v u yields 

\t+ I l w «l 2+ I QUe] k v t uk\ 2 = I iCttVu + [ Kvtt+H. (B.ll) 
1 at Jn Jn Jn Jr_ 

ii t 

We used the boundary condition (|B.4b[) in writing the boundary integral i. Thanks to 
Lemma I2TT1 we have that |u«|q = it- Hence, i = 7Z. Taking the time-integral of (|B.11[) and 
using that g > X e > 0, 

sup ||t>tt(i)||g+ / \\v t t\\j < [ K. (B.12) 
te[o,T] Jo Jo 

Step 2: Elliptic estimates for v t . A time-derivative of the equations (|B.4[) yields the 
following linear Neumann- type elliptic problem for i> t : 

~g[A t ]i([A t ] k r v uk ), = Q x in Q, (B.13a) 

gN j [Ae]i[A € ] s r v t , s = Ji on T, (B.13b) 

where the forcing functions Q\ and Ji are defined as 

Qx = jC t - v tt + {g[A t ]l) t ([A t ] k r v, k ) yJ +g[A t ]l(dt[A t ] k r v, k ) rJ , 

n = d t [W + d>{t)] - (gN j [A e ]i[A e Y r ).v, s . 
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Since e[I £ ]^j[[i e ]*9jt] is a uniformly elliptic operator, standard elliptic regularity theory 
provides that for s > 2 

INI* ^cdMl + w^wU + lMl-i), CB.14) 

with the constant C depending on the coefhcients of the elliptic and Neumann-type opera- 
tors. The fundamental theorem of calculus provides that 

lkllo<Ao + CT sup |M*)||<! < / K. 

te[o,T] Jo 

Thanks to the estimate (|B.12|) . we conclude that 



sup \\§ i (t)\\l + [ < [ k. 

;G[0,T1 Jo JO 



te[o,T] 

Since ji scales like Dv + Dv + Vt, the fundamental theorem of calculus provides that 

SUP |jl(*)|g. B + f |jl|?. 5 < f Tl+ SUp \\v t (t)\\l 
t€[0,T] Jo Jo te[0,T] 

Since interpolation and Young's inequality provide that \\v t \\l < Call^tHo + ^11^*111) we con " 
elude by use of the fundamental theorem of calculus that 

sup / Ulll.6 < / K. 

te[o,T] Jo Jo 

It therefore follows from the inequality (|B. 14|) that 

sup \\v t (t)\\l + [ \\v t \\l< [ It. (B.15) 
te[o,T] Jo Jo 

Step 3: Concluding the proof of Lemma IB. 21 By repeating Step 2, we infer the 
following elliptic estimate for v: 

sup \\v(t)\\t+ [ \\v\\l< [ 1i. (B.16) 



te[o,T] 

The sum of the inequalities (|B.12|) . (|B.15I) and (|B.16|) establishes the inequality (|B.10|) . □ 

Lemma B.3. Let v be the solution of (|B.4[) determined by v £ Ct{M). The solutions map 

&:C T (M)^C T {M): v^v 

is a well-defined map for some T — T K (e/j,) and possesses a unique fixed-point. 

Proof. We recall that v € is unique and satisfies |(X1)| by Lemma IB. 21 . Setting M. = 
JVo + 1, we take 5 so that i5||S||x r < \ an d let T — T K (e/i) be sufficiently small so that the 
inequality given in Lemma IB.2I reads 



\Hk T <M. 

Hence, for some T = T K (e^t) > 0, the solutions map & is well-defined. 

Letting vi € Ct{M.) for I = 1, 2, we set fji = e + f Q vi in f2 and f] e i = e + J Q v e i on T. We 
define f e j to be the solution of the following time-dependent elliptic Dirichlet problem: 

AQ = Afji in O, (B.17a) 

Cei = fj el on T. (B.17b) 
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We define the following e-approximate Lagrangian variables: 

A e i = [DC e ;]"\ Ji = dctD( el , a e i = J e iA e i, [g e i} a p = Ceha <el>/9 , VSTind = [a e i] T N. 
For I = 1, 2, we also take the quantities g~i, Ki, to be respectively formed via the definitions 



(|B.5I) , (|B.6[) , (|B.7|) with v — vi . Letting vi denote the solution of (|B.4|) formed using v — vi , 
we have that the difference 

w = V\ — V 2 

satisfies 

w t -Q 1 [A el }i([A el ]*w, k ), j =F inOx(0,T K ( eM )], (B.18a) 
g 1 N^[A £l Y r [A el ^w lS = G on V x (0,T K (e M )], (B.18b) 
w\ t =o = on Q, (B.18c) 
where the vector field F appearing in the right-hand side of (IB.18a[) is given by 

F = Ki - K 2 + gMeiYr([Aei} k r v 2 ^)^ - 92 [A e2 }i ( MS ,fc ) ,j , (B.19) 
s . ' 

3 

and the vector field G appearing in the right-hand side of (|B.18b[) is given by 

G = ht-h 2 l - [QiNi[A el ]i[A el ]*v 2 ,k +Q2N^A e2 y r [A e2 } k r v 2 , k ] . (B.20) 

" V ' 

j 

Testing two time-derivatives of (|B.18ap against ro (t in the L 2 (Sl)-inner product and inte- 
grating by parts in the integrals - f Q d^(gt[A e i]^.(JA e i]^w,k ),j )w u and J n dttw tt yields 

' " r \w tt \ 2 + [ d?(Q 1 [A el } J r [A el ]*w, k )wu, j + f d^({gi[A tl y r ), J [A tl }';w,k)wu 



2dt 



- - h%)ttw t t 



2 
't 

n 



Of 



w tt 



Ei -Ki - - (Q2{A e2 y r ), 3 [A t2 t)v 2lk 

- [ d?[(g 1 [A el y r [A el ] k r - g 2 [A e2 y r [A e2 ] k r )v 2 , k ]w tt , J ■ (B.21) 
Jn 

We have used that two time-derivatives of the boundary condition (|B.18b[) is equivalent to 

df^N^A^yiA^w^ +j] = - hi) tt . 

Since 9 2 (gi — g 2 ) — J * df (g~i — g 2 ) scales like J Q T D(v\ — v 2 ) t t, we infer that the time-integral 
of (jB~2T|) yields 

sup \\w tt (t)\\ 2 + [ \\wuWl < 5 [ ||«i-«2||x T +^ T ll«i-^||x T - (B.22) 
*e[o,T] Jo Jo 

We use the estimate (|B.22I) and follow Steps 2 and 3 of the proof of Lemma IB. 21 to obtain 
estimates for w t in L 2 (0, T; H 3 (£l)) and w in L 2 (0, T; H 5 (£l)). Thus, 

WL<<5 / T ||wi-«2|| Xr +C^^II«i-«2llx r - (B.23) 
Jo 

We recall that iu = vi — v 2 . The inequality (|B.23|) therefore provides that 6 is a contraction 
mapping for sufficiently small S > and T — T K (efi) > 0. Hence, the solutions mapping 6 
possess a unique fixed-point v € Ct(A4). □ 
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B.4. The proof of Lemma lB.ll The proof of Lemma lB.ll is complete since the fixed-point 
v € Ct{M.) of Lemma [B . 3 1 verifies | (XI ) | and is the unique solution of the ^(-problem ()5.52|) . 

B.5. The proof of Proposition [B7T1 Wc infer from the proof of Lemma IB.ll that con- 
sidering two more time-derivatives in the functional framework of our fixed-point scheme 
establishes Proposition IB.ll 



Appendix C. Equivalence of the ace-problem (I5.6[) and the heat-type 

K6-PROBLEM (|5.49l) 

In this appendix, we prove the following 

Lemma C.l. The ne-problem (|5.6j) and the heat-type ne-problem (|5.49|) are equivalent. 

Proof. Sections 15.3.11 and 15.3.21 provide that a solution of the «;e-problem (|5.6p satisfies the 
heat-type Ke-problem (|5.49[) . We now establish the converse. Using the identity (I5.29P stat- 
ing — [A e ]l [[A e ] k v,k ],j — cur l( £ CUY ^ t v ~ [Ae] s (div^ v), s , the nonlinear heat-type equations 
(|5.49a[) are equivalently written as 

v t + gcurl^ curl^ v — g[A e ] s (div^ v), s = K. 

The identity 

-g[A t ] s (dW L v), s = -[A e ] k (gdiv L v), k + div^ v[A t ] k g, k , 
and the definition (|5 . 32[) of JC imply that the equations (|5.49a[) are equivalently written as 

v t - [A e ] k (gdiv; v - 2p J~ 1 ), k = -pcurb [curb v - (curl u Q + s.ji / d t [A e ]jV l lS )]. 

Jo 

(C.l) 

Using the identity curl^ v = curl no + f Q 9t(curb v) = curluo + f 9t(e.ji[A e \jV' 1 , s ), 

/ curl^ v t = cur 1^ v - (curl u + e.ji / d t [A e ] s j v' l , s ). (C.2) 
Jo Jo 

The identity (|C2[) implies that the equations (|C.1[) are equivalently written as 

v t - [A e } k (gdiv; v-2p J- 1 ) lk = -gcuvh [ [ curb v t ]. (C.3) 

Jo 

Applying the e-approximate Lagrangian curl operator curl^ to (IC.3j) yields 

curb v t + curb (gcurb / curb v t ) — in Q. (C4) 
6 Jo 

We recall c[t) defined in (|5.50[) . Using the definition (|5.46l) of h, we equivalently have that 

2 .a 

c(t) = V ^d?lgN3[A e y r [A e y r v,s -g[bcnn + b d 



a'. 

a=0 

2 t a 

- ^ — t dt[h cm \ + h div + g^[v, a/3 ■n e ]n e ]|t=o- 
*■ — ' a! 

a=0 
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By the identities (|5.39|) and (|5.48p . we have that 
2 t a 

c(t) = 22 —d^[gy/f e J £ x (curl c - e v) x h e + g^/fj t ^div^ v)n e ]\ t=0 

a=0 °" s ' 

eN3[A e y r [A e ]sv, 3 -gb 

- -r^curi + gJ ( ~ 1 v / ir(div c - e v)n]\ t=0 

a=0 ^ ^ ' 



2 r 



= ^ ^r9 t [eV^J f e ~ 1 (curb w) x n e - ft cur l]|t=o- 
' — ' a! s 

a=0 

Recalling that h CUT \ = g^/glJ^ 1 (curlu + e.ji f Q dt[A e ]jV i , s ) x n e , we conclude that 

2 ft 

c(t) = V — [pv/^J'^curb e - curluo - e.^ / d t [A e ]* i v i » ) x n £ ]| t =o- 
By the identity (|C.2j) for J Q curl,; i)t, we have established that 

c(t) = ^ ^c^VW^ /*curl 6 «t) x n e ]| t=0 . (C.5) 

a=0 °' ^° 

We will now verify that c(t) = 0. With (J* curl^ Vt)\t=o — 0, we equivalently write (|C.5|) as 



c(t) = tKp ^g e curl vi x AT + y ftp V <?e [curl v 2 + e.- )i a t [A e ]|(0)v l 1 , s ] x N. 

According to (I5.10[) , we have that vi = D(np div m — 2po)- Hence, curlvi = 0. Using the 
identity d t [A e ]f\ t=0 = [A e ]^v r , s )\ t=0 = we conclude that 

t 2 I— 

C W = y«PoVffe[curlv 2 - e.jiUQ,j v\, s ] x AT. 
According to (|5.10p . v 2 = Dd t {np J t - 2p J e _1 )| t=0 + d t [A e ] k \ t=0 d k (Kp div u - 2p ). Thus, 
curlv 2 = -e.ji [u$,i (np divu - 2p ),fc ] ,j = -e. jl u k )n (npo divu - 2p ), jk = -e.^u^j vj, fe . 
That is, 

curlv 2 = E.jjWoy v i>« > 
so that the identity e(t) = ^-KpoVffTIcurl v 2 — e.jiV,Q,j v\, s ] x A" implies that 

c{t) = 0. 

The boundary condition ((5.49b)) is therefore equivalently written as 

QN j [A e ]i[A e ] s r v, s = 



hurt + Qb + — plJr 2 ~ &(*) + vgfCoaP 'fie + SS P i )a p -fl t fl e . (C.6) 
Po L J 

Taking the scalar product of (|C.6|) with g~ 1 8( e yields 

v /^j7 1 [(curlf e v) x h e ] ■ d( e = i/^J e _1 [(curlno + e.ji / d t [A e ]jV\ s ) X n e ] ■ 9(e ■ 



o 
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Using that [(curl^ v) x n e ] ■ n e = and /i C uri • n £ — 0, we have that 

(curb v) x n £ = (curlw +£-ji / 9t[-A e ]jU l , s ) x ri e on T. (C7) 

Jo 

We note that (IC.7|) is equivalently stated as (J Q curb u 4 ) x n e = 0. Hence, by the equation 
(|C.4[) , we record that J* curl,; v t satisfies 

curl^ Vt + curh^ (pcurl^ / curl^ fit) = in 51 x [0, T], 
' •/() 



(C.8a) 



( / curl fe v t ) x n c = on b x [0, T], (C.8b) 



( / curb v t )\t=o = on 17. 
Jo 



(C.8c) 



Testing (|C.8a[) against J Q curl^ w t in the L 2 (f2)-inner product and integrating by parts with 
respect to the operator curl^ in J a curl,; (pcurl^ f Q curl^ w t ) J Q curl^ v t yields 

1 d 



2 dt 



|| / curl^ u t ||o + || curb / curl^w t | 
Jo Jo 

+ I A " ,1, ; ,ji {ovnv\^ I curl c - e ^) fc (/ curl^ v t y = 0. (C.9) 



With the boundary condition (|C.8b|) . J r N s [A e ]j£ijk(gcuil^ f Q curb fit) k (f curl^ v t ) 1 = 

J r y/TjlJ^ 1 [( J Q curl^ U{) x n] ■ curl^ f Q curl^ v t ) = 0. Integrating the identity (|C.9I) in time 

from to t G (0, T] produces || f Q curl^ j5 4 ||q = 0. So according to (|C.8a|) . || curl^ # t || = 0. 
The regularity of v stated in Proposition 15 .41 then provides that almost everywhere in fi, 



curl^ Vt = 0. 

The equations (IC3|) are thus written as 

v t - [A e ] k (K Po J e div £c v - 2 PQ J~ l ) lk = 0. (C.lOa) 

As the solution v of the heat-type Ke-problem (|5.49l) verifies the e-approximate bagrangian 
vorticity equation (j5.15[) , we infer that the boundary condition (|C.6|) of the heat- type prob- 
lem is equivalently written as 

K P 2 div ie v = p 2 J- 2 - e (t) + agfCe, a p -n e + K~gfv, aP -n e . (C.lOb) 

The equations (IC.10[) are equivalent to the momentum equations and boundary condition 
of the /te-problem (|5.6[) . □ 
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List of Notation 

u Eulerian velocity [2] 

p Eulerian pressure ?? 

p Eulerian density [2] 

u Q Initial velocity ?? 

po Initial density ?? 

a Surface tension ?? 

(3 A parameter to the equation of state p = ap 1 — j3 for 7 > 1 [2] 

H (77) Twice the mean curvature of the moving surface 77 (r) ?? 

rj The particle flow map [3] 

v Lagrangian velocity [3] 

/ Lagrangian density [3] 

A Inverse of the deformation tensor Dr\ [3] 

J Jacobian determinant of the deformation tensor Drj [3] 

a Cofactor matrix of the deformation tensor Drj [3J 

e The identity function e(x) = x [4] 

D The three-dimensional gradient operator [7] 

8 The surface gradient operator [7] 

divjjjCurl,; The Lagrangian divergence and curl operators [7] 

n The outward-pointing unit normal to the moving surface r)(T) 

g The surface metric induced by the moving surface rj (L) [9] 

go The surface metric of the initial surface L M 

^fg The square root of the determinant of the metric g [9] 

A g The Laplace-Beltrami operator [9] 

|| • || s The norm of the Hilbert space H S (Q) M 

I • \ s The norm of the Hilbert space H S {T) [9] 

N The outward-pointing unit normal to L [10] 

v A solution of the K-problem (|3.1[) [15] 

?y The Lagrangian map of the solution v to the K-problem [15] 

E K (t) The higher-order energy function for solutions v to the re-problem . . [15] 

V, Ao, 1Z The generic system of constants used in the k- independent estimates [15] 

A c The horizontal-convolution operator for e > [27] 

v A solution of the ree-problem ()5.6|) [28] 

E e (t) The higher-order energy function for solutions v to the ree-problem . [28] 

V,Nq, 1Z The generic system of constants used in the e- independent estimates [25] 

v A solution of the /^-problem (|5.52[) [35] 

The higher-order energy function for solutions v to the /^-problem __. 

(I5321 ^ 

VjAfo, 1Z The generic system of constants used in the p- independent estimates [35] 

v A solution of the cr-problem ([7T2| [42] 

The higher-order energy function for solutions v to the surface tension 

problem 

P, Ao, 7?. The generic system of constants used in the cr-independent estimates [42] 

Xt The Hilbert space used in the fixed-point scheme [53J 

Ct(M) A closed, bounded, convex subset of Xt [63] 

v An arbitrary vector in Ct(-M) [64] 

V, Ao, Ti- The generic system of constants used in the fixed-point scheme .... [65l 
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